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25.1. Introduction 

We have already discussed in chapters 7 and 8, the 
graphical methods to determine velocity and acceleration 
analysis of a mechanism. It may be noted that graphical 
method is only suitable for determining the velocity and 
acceleration of the links in a mechanism for a single position 
of the crank. In order to determine the velocity and 
acceleration of the links in a mechanism for different 
positions of the crank, we have to draw the velocity and 
acceleration diagrams for each position of the crank which 
is inconvenient. In this chapter, we shall discuss the analytical 
expressions for the displacement, velocity and acceleration 
in terms of general parameters of a mechanism and 
calculations may be performed either by a desk calculator 
or digital computer. 


19. Analytical Synthesis for 
slider Crank Mechanism 
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25.2. Computer Aided Analysis for Four Bar Mechanism (Freudenstein’s 
Equation) 

Consider a four bar mechanism ABCD, as shown in Fig. 25.1 (a), in which AB = a, BC = b, 
CD = c, and DA = d. The link AD is fixed and lies along X-axis. Let the links AB (input link), BC 
(coupler) and DC (output link) make angles 0, p and <|) respectively along the X-axis or fixed link 
AD. 




( b ) Components along X-axis and T-axis. 

Fig. 25.1 


The relation between the angles and link lengths may be developed by considering the 
links as vectors. The expressions for displacement, velocity and acceleration analysis are derived 
as discussed below : 


1. Displacement analysis 

For equilibrium of the mechanism, the sum of the components along X-axis and along 
T-axis must be equal to zero. First of all, taking the sum of the components along X-axis as shown 
in Fig. 25.1 ( b ), we have 

flcos9 + fecosP-ccos(j)-c/= 0 ...(/) 

or ZjcosP = ccos(]) + d-acos0 

Squaring both sides 

b~ cos 2 P = (c cos (]t + d — a cos 0) 2 

- c 2 cos 2 (p + d 2 + 2c d cos (]) + a 2 cos 2 0 

-2a ccos<j)cos9-2a dcosQ . . . (ii) 

Now taking the sum of the components along T-axis, we have 

a sin 9 + b sin P - c sin <j) = 0 ... (ttt) 

or b sin P = c sin (]) - a sin 0 

Squaring both sides, 

b 2 sin 2 P = (c sin (f) — a sin 9) 2 

= c 2 sin 2 (]) + a 2 sin 2 9 — 2a csin(|)sin9 . . . (iv) 

Adding equations (it) and (tv), 

b 2 ( cos 2 P+sin 2 P) =c 2 (cos 2 + sin 2 §) + d 2 +2 c d cos(]) + a 2 (cos 2 0 + sin 2 9) 




-2 a c(cos<j)cos0 + sin(()sin0)-2 a d cos0 
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b l = c z +d l +2 c d cost]) + a -2 a c(cos(])cos9 + sin(|)sin0)-2 a d cos0 

2 a c(cos(])cos0 + sin(])sin0) = a 2 -b 2 +c 2 +d 2 +2 c d cos(])- 2a d cos0 

■ „ . , . „ a 2 -b 2 +c 2 +d 2 d , d a 

cos (pcos 0 + sin <p sin 0 =-1— cos <p-cos 0 

2 a c a c 


Let 


— = k \; — = k 2 ; and 


a -b +c +d 
2 a c 


— 


■ (v) 

(vi) 


Equation (v) may be written as 

cost)) cos 0 + sin (j> sin 0 = Aq cos t() — A: 2 cos 0 + A: 3 ... (vii) 

or cos (<|) - 0) or cos(9-()>) = Aq cost))- k 2 cos9 + k 3 

The equation (vii) is known as Freudenstein’s equation. 

Since it is very difficult to determine the value of <|) for the given value of 0, from 
equation (vii), therefore it is necessary to simplify this equation. 

From trigonometrical ratios, we know that 

. 2 , 


2tan((j>/ 2) 
1+ tan 2 (<j)/2) 


and cos tb = 


1 — tan (()> / 2) 
1 + tan 2 (())/ 2) 


Substituting these values of sintj) and cost]) in equation (vii), 

1 — tan 2 (d) / 2) „ 2tan(tb/2) . „ 

- — 1 XCOS0 + -- 1 - xsin 0 


1 + tan (tf> / 2) 


1 + tan (tf) / 2) 


_, l-tan 2 (t|)/2) 

— x-—- ^2 cos 0 + 


1 + tan (tf) / 2) 
cos0 [1 - tan 2 (t|)/ 2)] + 2sin 9 tan (tf) /2) 


= k, [I - tan 2 (t|) / 2)] - k 2 cos 0 [1 + tan 2 «)/ 2)] + k 3 [1 + tan 2 (t^/ 2)] 


cos 0-cos 0 tan 2 (tf)/ 2) + 2sin 9 tan (tf>/ 2) 


= k x -k x tan 2 (t])/ 2 )-k 2 cos0-k 2 cos0tan 2 (t|)/2) + k 3 +k 3 tan 2 (tf>/2) 
Rearranging this equation, 

—cos0tan 2 (t|)/2) + /: 1 tan 2 (t|)/2) + k 2 cos9tan 2 (t|)/2) — k 3 tan 2 (t|)/2) + 2sin0tan (tf)/2) 

= —cos0 + /q — k 2 cos0 + k 3 

— tan 2 (tf)/2) [cosO-kj —k 2 cos0 + k 3 ] + 2sin0tan (t|)/2)-k 1 -k 3 +cos0(l + k 2 ) = 0 

[ (l-k 2 )cos0 + k 3 -k\ J tan 2 t|)/2 + (-2sin 0)tant|)/2 + [k 1 +k 3 -(l + k 2 ) cos 9] =0 

(By changing the sign) 

Atan 2 ((J)/2) + i? tan((|) / 2) + C = 0 ... (viii) 



or 
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where A = (l — ) cos 0 + ^3 +k\ 

B = — 2 sin 0, and 

C = /q + k ^ — (1 + ^ 2 ) cos 0 


■ ■ ■ (ix) 



Inner view of an aircraft engine. 

Note :This picture is given as additional information and is not a direct example of the current chapter. 


The equation (vtti) is a quadratic equation in tan(<j)/2). Its two roots are 


tan ((f)/ 2 ) = 


-B±\Ib 2 -4AC 


2A 


or 


(f) = 2 tan 1 


-B±^B 2 -4AC 
2A 


■ ■ ■ to 


From this equation (x), we can find the position of output link CD (i.e. angle (f>) if the 
length of the links (i.e. a, h. c and d) and position of the input link AB (i.e. angle 0 ) is known. 

If the relation between the position of input link AB (i.e. angle 0) and the position of 
coupler link BC (i.e. angle (3 ) is required, then eliminate angle (|) from the equations (i) and (Hi). 
The equation ( 1 ) may be written as 

ccos(|) = acos0 + bcos(3 — d ... (xi) 

Squaring both sides, 

c 2 cos 2 (|) = a 2 cos 2 Q + b 2 cos 2 (3 + 2 a bcos0cos(3 


+ d 2 -2 a dcos 0-2 b d cos(3 
Now equation (iii) may be written as 
c sin (|) = a sin 0 + b sin (3 


. . . (xii) 



. . . (xiii) 
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Squaring both sides, 

c 2 sin 2 (f) = a 2 sin 2 Q + b 2 sin 2 (3 + 2 a b sin 0 sin(3 ... (xiv) 

Adding equations (xii) and (xiv), 

c 2 ( cos 2 (|) + sin 2 0 ) = a 2 (cos 2 0 + sin 2 0)+b 2 (cos 2 (3 + sin 2 (3) 

+ 2aZ?(cos0cos|3 + sin0sin(3) + (i 2 -2 ad cos 0-2 b d cos (3 


c = a +d +2 a i?(cos0cos(3 + sin0sin(3) 

+d 2 -2adcosQ-2b d cos(3 

2a/?(cos0cos(3 + sin0sin(3) = c 2 -a 2 -b 2 -d 2 +2a d cos0 + 2 b d cos(3 

„ D ■ « • o c 2 -a 2 -b 2 -d 2 d . d _ 

>s0cos[3 + sin0sin(3 =- 1 —cos0H—cos|3 


Let 


2ab 


d , d , . c — a —o — a , 

- = ki;- = k 4 ; and -—-= ^ 5 

a b 2a b 


Equation (jcvi) may be written as 

cos 0 cos (3 + sin 0 sin (3 = &q cos (3 + k 4 cos 0 + k 5 
From trigonometrical ratios, we know that 


2tan((3/2) 
1+ tan 2 ((3/2)' 


and cos (3 = 


1 — tan 2 ([3 / 2 ) 
1 +tan 2 ((3/2) 


Substituting these values of sin (3 and cos (3 in equation (xvii), 
cos 0 


1 -tan 2 ((3/ 2 ) 

+ sin 0 

2 tan( (3 / 2) 

1 + tan 2 ((3 / 2 ) 

1 + tan 2 ((3 / 2 ) _ 



= k. 


1 — tan 2 ([3 / 2 ) 


+ k 4 cos 0 + ^5 


• . (XV) 

. (xvi) 

. (xvii) 


1 + tan 2 ((3 / 2 ) _ 
cos 0[1 — tan 2 ((3 / 2)] + 2 sin 0 tan( (3/2) 

= k\ |^ 1 - tan 2 ((3/ 2 )J + k 4 cos 0^1 + tan 2 ((3/2)J +k 5 ||l + tan 2 ((3/ 2 )J 

cos 0 — cos 0 tan 2 ((3 / 2) + 2 sin 0 tan( (3/2) 

= /q -/q tan 2 ((3/2) + A : 4 cos0 + fc 4 cos 0 tan 2 ((3/2) 

+k 5 +k 5 tan 2 ((3/2) 

-cos0tan 2 ((3/2) + /q tan 2 ((3/2) — A : 4 cos0tan 2 ((3/2 )-/: 5 tan 2 ((3/2) 

+ 2sin0tan((3/2)-/q -k 4 cos Q-k 5 +cos0 = O 
— tan 2 ((3/2)[(At 4 + 1)cos0 + £ 5 -/q] + 2sin0tan(|3/2 )-[(/: 4 -l)cos0 + /: 5 + /q] =0 

















Chapter 25 : Computer Aided Analysis and Synthesis of Mechanisms • 1007 


[(k 4 +l)cos 0 + fc 5 -^ 1 ]tan 2 (p/ 2 ) + (-2sin0)tan((3/2) + [(^ 4 -l)cos 0 + /: 5 + /q] =0 

(By changing the sign) 


where 


Dtan 2 (p/2) + £tan(p/2) + F = 0 

D = (k 4 +1) cos 0 + ^5 —k\, 

E = —2 sin 0, and 
F = [(k 4 - 1)cos0 + k$ + /q] 

The equation (xviii) is a quadratic equation in tan((3/2). Its two roots are 
-E±^E 2 -4DF 


(xviii) 


.. ( xix ) 


tan((3/2) = 


(3 = 2 tan 


2D 

-E±\Ie 2 -4DF 


2D 


(xx) 


From this equation (xx), we can find the position of coupler link BC (i.e. angle P ). 

Note: The angle a may be obtained directly from equation (i) or (hi) after determining the angle ()). 

2. Velocity analysis 

Let COj = Angular velocity of the link AB = dQ/ dt , 

tt >2 = Angular velocity of the link BC = c/p / dt, and 

CO 3 = Angular velocity of the link CD = c/<|)/ dt. 

Differentiating equation (i) with respect to time, 

c /0 „ c/p . c/(j) 

— a sin 0 x- b sin P x — + c sin <j) x — = 0 

dt dt dt 

or -aCOj sin 0 —/?(02 sinP+cCOj sin(|) = 0 

Again, differentiating equation (Hi) with respect to time, 

„ C/0 , „ £/P , C/(J) 

a cos 0 x-h b cos P x ——— c cos <() x —— = 0 

dt dt dt 

or aco 1 cos0 + / 7 CO 2 cos P -c ® 3 cos( l ) = O ... (xxii) 

Multiplying the equation (xxi) by cos P and equation (xxii) by sin P, 

—ctCOj sin0cosP-bt 02 sinPcosP + c (03 sin(|)cosP = 0 ... (xxiii) 

and aCOj cos0sinp + (0 2 cosPsinp — cCO 3 cost^sinP = 0 ... (xxiv) 

Adding equations (xxiii) and (xxiv), 

a (Oj sin(P - 0) + c co 3 sin (<|> - P) = 0 


(xxi) 
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Again, multiplying the equation (xxi) by cos<j> and equation ( xxii ) by sin(|), 

-aCOj sin0cos(j)-bo>2 sin (3 cost])+ cco 3 sin(j)cost]) = 0 ... (xxvi) 

and a COj cos 0 sin cf) + b co 2 cos P sin (]) - c (O3 cos (]) sin (]) = 0 ... (xxvii) 

Adding equations (xxvi) and (xxvii), 

a (Oj sin (tj) - 0) + b co 2 sin(4> — P) = 0 
-a sin( (|) — 0) 


C0 2 = 


... (xxviii) 


Z? sin(<|) — P) 

From equations (xxv) and (xxviii), we can find ©3 and co 2 , if a, b, c, 0, <]). p and © ( are 

known. 

3. Acceleration analysis 

Let a! = Angular acceleration of the link AB = d®\ / dt, 

a 2 = Angular acceleration of the link BC = c/© 2 / dt , and 
a 3 = Angular acceleration of the link CD = d C0 3 / dt. 

Differentiating equation (xxi) with respect to time. 


J0 d® x 

©! cos 0x-h sin 0x —- 

dt dt 


— b 


0 ^P . 0 d® 2 

C0 9 cosPx —+ sinPx—^ 
dt dt 


. tM> . , d®T. 

©3 cost]) x — +sint]) x—— 
dt dt 


= 0 


—a cof cos0-asin0a! —box) cosP-Z)sinPa 2 
+CCO3 cos(]) + c sin(]>a 3 = 0 

Again, differentiating equation (xxii) with respect to time. 


dQ d©! 

©! x - sin 0 x-f cos 0 x —- 

dt dt 


+ b 


d , dv du 

-( UV) = UX -h U X 


. „ ^P 0 d ® 2 

© 7 x - sin Px — +cosPx—- 
dt dt 


dx 


dx 


(xxix) 


. , dti) , d® 3 
©3 x - sin t]) x — + cos ({) x — 
dt dt 


= 0 


- a©! sin 0 + a cos 0a t — b © 2 sin P + b cos P a 2 
+ca^ sint])-ccos(])a3 =0 

Multiplying equation (xxix) by cost]), and equation (xxx) by sint]), 
—a ©j 2 cos 0cos § — a a! sin 0cost])- b®\ cosPcost]) 


(xxx) 


—Z>a 2 sin Pcost]) + c ©3 cos^ t]) + ca 3 sin tpcos t|) = 0 



... (xxxi) 

















n n 
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and 


and 


—aCOj 2 sin 0 sin (]) + a oq cos9sin(|)— bco 2 sin(3sinc() 

+/?a 2 cos(3sin(|) + cC03 sin 2 (|)-ca3 cos(|)sin(t) = 0 (xxxii) 

Adding equations (xxxi) and (xxxii), 

—a (0^ (cos (]> cos 9 + sin (]) sin 0) + a oq (sin (]) cos 0 — cos (]) sin 9) 

—b ($2 (cos (]) cos P + sin (]) sin (3) + b a 2 (sin (|> cos P - cos (f>sin P) 

+c CO 2 (cos 2 (|) + sin 2 tf)) = 0 

-a CO 2 cos((|)-0) + fla 1 sin (cp — 0) — &co| cos((])-P) + /;>a 2 sin ((|) — P) + c-cof =0 


a 2 = 


-aa x sin (4 - 0) + a (of cos((])-0) + bco 2 cos((|)-P)-c(03 


Z?sin (<|> —P) 

Again multiplying equation (xxix) by cos P and equation (xxx) by sin P , 
-act)! 2 cos0cosp-aoq sin0cosP- b to 2 cos 2 P~ba 2 sinpcosp 

+c CO 2 cos (|) cos P + c a 3 sin (|) cos P = 0 
-aCOj 2 sin0sinP + aoq cos0sinp-ba> 2 sin 2 P+ba 2 cosPsin P 


+c(0 2 sin(()sinP-ca 3 cos(])sinP = 0 


Adding equations (xxxiv) and ( xxxv ), 


(xxxiii) 


(xxxiv) 


... (xxxv) 


—a to 2 (cos P cos 9 + sin P sin 0) + a oq (sin P cos 0 — cos P sin 0) — b (0 2 (cos 2 P + sin 2 P) 
+ c CD 2 (cos (|)cos P + sin (|) sin P) + c a 3 (sin (])cos P - cos (psin p) = 0 
-a to 2 cos (P - 0) + a sin (p - 0) -i> co^ + c to 2 cos ((|)- P) + c a 3 sin ((|) - P) = 0 


a 3 ~ 


-aa { sin(p-0) + aco 2 cos(P-0) + bco 2 -c CO 2 cos(tf> — P) 
csin ((p — p) 


... (xxxvi) 


From equations (xxxiii) and (xxxvi), the angular acceleration of the links BC and CD (i.e. 
a 2 and a 3 ) may be determined. 


25.3. Programme for Four Bar Mechanism 

The following is a programme in Fortran for determining the velocity and acceleration of 
the links in a four bar mechanism for different position of the crank. 

PROGRAM TO FIND THE VELOCITY AND ACCELERATION IN A FOUR-BAR 
MECHANISM 

DIMENSION PH (2), PHI (2), PP (2), BET (2), BT (2), VELC (2), VELB (2), ACCC (2), 

ACCB (2), Cl (2), C2 (2), C3 (2), C4 (2), B1 (2), B2 (2), B3 (2), B4 (2) 

READ (*, *) A, B, C, D, VELA, ACCA, THETA 
PI = 4.0 * ATAN (1.0) 

THET = 0 
IHT = 180/THETA 
DTHET = PI/IHT 
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DO 10 J = 1, 2 * IHT 

THET = (J - 1) * DTHET 

AK = (A*A-B*B+C*C + D*D)* 0.5) 

TH = THET * 180/PI 

AA = AK - A * (D - C) * COS (THET) - (C * D) 

BB = - 2.0 * A* C * SIN (THET) 

CC = AK - A * CD + C) * COS (THET) + (C * D) 

AB = BB * * 2 - 4 * AA * CC 
IF (AB . LT . 0) GO TO 10 
PHH = SQRT (AB) 

PH (1) = - BB + PHH 
PH (2) = - BB - PHH 
DO 9 I = 1, 2 

PHI (I) = ATAN (PH (I) * 0.5/AA) * 2 
PP (I) = PHI (I) * 180/PI 

BET (I) = ASIN ((C * SIN (PHI (I)) - A * SIN (THET)) / B) 

BT (I) = BET (I) * 180/PI 

VELC (I) = A * VELA * SIN (BET (I) - THET) / (C * SIN (BET (I) - PHI (I))) 

VELB (I) = (A * VELA * SIN (PHI (I) - THET) ) / (B * SIN (BET (I) - PHI (I)))) 

Cl (I) = A * ACCA * SIN (BET (I) - THET) 

C2 (I) = A * VELA * * 2 * COS (BET (I) - THET) + B * VELB (I) * * 2 
C3 (I) = C * VELC (I) * * 2 * COS (PHI (I) - BET (I) ) 

C4 (I) = C * SIN (BET (I) - PHI (I)) 

ACCC (I) = (Cl (I) - C2 (I) + C3 (I) ) / C4 (I) 

B1 (I) = A* ACCA* SIN (PHI (I) - THET ) 

B2 (I) = A * VELA * * 2 * COS (PHI (I) - THET ) 

B3 (I) = B * VELB (I) * * 2 * COS (PHI (I) - BET (I) ) - C * VELC (I) * * 2 
B4 (I) = B * (SIN (BET (I) - PHI (I)))) 

9 ACCB (I) = (B1 (I) - B2 (I) - B3 (I)) / B4 (I) 

IF (J . NE . 1) GO TO 8 

WRITE (*, 7) 

7 FORMAT (4X,’ THET’, 4X,’ PHI’, 4X,’ BETA’, 4X,’ VELC’, 4X,’ VELB’, 4X,’ ACCC’, 4X,’ 
ACCB’) 

8 WRITE (*, 6) TH, PP (1), BT (1), VELC (1), VELB (1), ACCC (1), ACCB (1) 

6 FORMAT (8F8 . 2) 

WRITE (*, 5) PP (2), BT (2), VELC (2), VELB (2), ACCC (2), ACCB (2) 

5 FORMAT (8X, 8F8 . 2) 

10 CONTINUE 
STOP 
END 

The various input variables are 

A, B, C, D - Lengths of the links AB, BC, CD, and DA respectively in mm, 

THETA - Interval of the input angle in degrees, 

VELA = Angular Velocity of the input link AB in rad/s, and 
ACCA = Angular acceleration of the input link in rad/s 2 . 

The output variables are : 

THET = Angular displacement of the input link AB in degrees, 

PHI = Angular displacement of the output link DC in degrees, 

BETA = Angular displacement of the coupler link BC in degrees, 

VELC = Angular velocity of the output link DC in rad/s, 

VELB = Angular velocity of the coupler link BC in rad/s, 

ACCC = Angular acceleration of the output link DC in rad/s 2 , 

ACCB = Angular acceleration of the coupler link BC in rad/s 2 . 
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Example 25.1. ABCD is a four bar mechanism, with link AD fixed. The lengths of the links 


AB = 300 mm; BC = 360 mm; CD = 360 mm and AD = 600 mm. 

The crank AB has an angular velocity of 10 rad/s and an angular retardation of 30 rad/s 1 , 
both anticlockwise. Find the angular displacements, velocities and accelerations of the links BC 
and CD, for an interval of 30° of the crank AB. 

Solution. 


Given input : 


A = 300, B 

OUTPUT : 

= 360, C = 

360, D = 600, VA 

= 10, 

ACCA = 

-30, THETA = 

30 

THET 

PHI 

BETA 

VELC 

VELB 

ACCC 

ACCB 

.00 

- 114.62 

- 65.38 

10.00 

- 10.00 

- 61.67 

121.67 


114.62 

65.38 

10.00 

- 10.00 

121.67 

- 61.67 

30.00 

- 144.88 

- 82.70 

- 8.69 

- .84 

101.52 

181.43 


97.30 

35.12 

- .84 

- 8.69 

181.43 

101.52 

60.00 

- 166.19 

- 73.81 

- 6.02 

6.02 

38.02 

77.45 


106.19 

13.81 

6.02 

- 6.02 

77.45 

38.02 

90.00 

174.73 

- 47.86 

- 8.26 

12.26 

- 180.18 

216.18 


132.14 

- 5.27 

12.26 

- 8.26 

216.18 

- 180.18 

270.00 

- 132.14 

5.27 

12.26 

- 8.26 

- 289.73 

229.73 


- 174.73 

47.86 

- 8.26 

12.26 

229.73 

- 289.73 

300.00 

- 106.19 

- 13.81 

6.02 

- 6.02 

- 113.57 

- 1.90 


166.19 

73.81 

- 6.02 

6.02 

- 1.90 

- 113.57 

330.00 

- 97.30 

- 35.12 

- .84 

- 8.69 

- 170.39 

- 49.36 


144.88 

82.70 

- 8.69 

- .84 

- 49.36 

- 176.39 


25.4. Computer Aided Analysis For 

A slider crank mechanism is shown in Fi 
rod BC of length b. Let the crank AB of radius 



Fig. 25.2 Slider 


Slider Crank Mechanism 

;. 25.2 (a). The slider is attached to the connecting 
3 rotates in anticlockwise direction with uniform 



(■ b ) 

mechanism. 


angular velocity C0j rad/s and an angular acceleration a.| rad/s 2 . Let the crank makes an angle 0 
with the X-axis and the slider reciprocates along a path parallel to the X-axis, i.e. at an eccentricity 
CD = e, as shown in Fig. 25.2 (a). 
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The expressions for displacement, velocity and acceleration analysis are derived as 
discussed below : 

1. Displacement analysis 

For equilibrium of the mechanism, the sum of the components along X-axis and along Y- 
axis must be equal to zero. First of all, taking the sum of the components along X-axis, as shown in 
Fig. 25.2 (ft), we have 

a cos 0 + ftcos(— (3) — x = 0 ... ((3 in clockwise direction from X-axis is taken - ve) 

or ftcosP = x — acosG ... (i) 

Squaring both sides, 

ft 2 cos 2 (3 =x 2 +a 2 cos 2 0-2xacos0 ... (ii) 

Now taking the sum of components along T-axis, we have 
b sin(- (3) + e + a sin 0 = 0 
or -ftsinP + e = asin0 

ft sin P = e-asin0 ... (iii) 

Squaring both sides, 

ft 2 sin 2 P = e 2 +a 2 sin 2 0 —2easin 0 ... (tv) 

Adding equations (ii) and (iv), 

ft 2 (cos 2 P+sin 2 P) =x 2 +e 2 +a 2 (cos 2 0 + sin 2 0)-2xacos0-2easin0 
ft 2 =x 2 +e 2 +a 2 -2xacos0-2easin0 
or x 2 +(—2acos0)x + a 2 -ft 2 +e 2 -2easin0 = 0 


or x 1 +k l x+k 2 =0 

where k x = -2acos0, and k 2 = a 2 -ft 2 +e 2 -2easin0 

The equation (v) is a quadratic equation in x. Its two roots are 

~k\ ±^k 2 —4k 2 

x = --- 

2 


... (v) 

... (VI) 


... (vii) 


From this expression, the output displacement x may be determined if the values of a, ft, e 
and 0 are known. The position of the connecting rod BC (i.e. angle P) is given by 


or 


sin (—P) = 


a sin 0 - e 
ft 

e - a sin 0 


P =sin 




e — a sin 0 


... (viii) 
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Note : When the slider lies on the X-axis, i.e. the line of stroke of the slider passes through the axis of 
rotation of the crank, then eccentricity, e = 0. In such a case, equations (vi) and (viii) may be written as 

k\ = -2flcos0 , and fc 9 =fl 2 -6 2 
_l i -asin0 


and 


i = sin 


2. Velocity analysis 

Let (Oj = Angular velocity of the crank AB = dQt dt, 

C0 2 = Angular velocity of the connecting rod BC = t/(3/ dt , and 

v$ = Linear velocity of the slider = dx/ dt. 

Differentiating equation (i) with respect to time, 

. „ d(3 dx . dQ 

6x-sinpx — =-ax-sin 0x— 


and 


dt dt 


dt 


dx 


-a ©[ sin 0 - & co 2 sin (3 —— = 0 
dt 

Again, differentiating equation (Hi) with respect to time, 
c/(3 a dQ 

0COS|3X —- = -ACOS0X- 

dt dt 

a ©! cos Q + b ©2 cos (3 = 0 

Multiplying equation ( ix ) by cos (3 and equation (x) by sin [3 , 

dx 

-flCOj sin0cos(3-6 co 2 sin(3cos(3—^xcos(3 = 0 

dt 

a tOj cos 0 sin (3 + b co 2 cos (3 sin (3 = 0 
Adding equations (xi) and (xii), 

dx 

a ©[ (sin (3 cos 0 - cos (3 sin 0)—^xcos(3 = 0 
dt 

dx 

a©! sin((3-0) = — xcos(3 
dt 

dx a©jsin((3-0) 
dt cos (3 

From this equation, the linear velocity of the slider ( p s ) may be determined. 


■ (ix) 

.. (x) 

■ (xi) 
(xii) 


(xiii) 


The angular velocity of the connecting rod BC (i.e. C0 2 ) ma Y be determined from equa¬ 
tion (x) and it is given by 
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3. Acceleration analysis 

Let a] = Angular acceleration of the crank AB = r/W] / dt , 

0-2 = Angular acceleration of the connecting rod = c/© 2 / dt, and 

a s = Linear acceleration of the slider = d~xl dt 2 
Differentiating equation (ix) with respect to time. 


n d 9 d C0i 

to, cos9x-hsin0x—- 

dt dt 


— b 


o d \3 c/© 2 

C0 2 cos p x — + sin p x 


dt 


d A x 


dt \ dt 2 


= 0 



The chain-belt at the bottom of a bulldozer provides powerful grip, spreads weight 
and force on the ground, and allows to exert high force on the objects to be moved. 
Note : This picture is given as additional information and is not a direct example of the current chapter. 


Differentiating equation (ur) with respect to time, 
+ b 


. . c/0 . do\ 

to, x-sm0x-i-cos0x—- 

dt dt 


<sh x—sinBx —+ cosBx- 

dt dt 


= 0 


a ^otj cos 0 - © 2 sin 0J + b |^a 2 cos (3 - ©2 sin (3J = 0 

Multiplying equation (xv) by cos (3 and equation (xvi) by sin (3 , 

—a j^cxj sin0cos(3 + co^ cos0cos pj — ||cx 2 sin P cos P + © 2 cos 2 pj 

- ^-^xcosp = 0 
dt 2 

aj^oq cos0sinp-©j 2 sin0sinpj + b|^a 2 cos psin P-©^ sin 2 pj = 0 


... (xvi) 


... (xvii) 


and 


... (xviii) 
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Adding equations (xvii) and (xviii), 

a j^OCj (sin (3 cos 0 — cos (3 sin 0) — (Of (cos (3 cos 0 + sin (3 sin 0) J 

-bco^(cos 2 p + sin 2 p)-^-^xcosp = 0 
dt 2 


i cq sin (P - 0) - a cof cos (P - 0) - b (£) 2 -— x cos p = 0 

dt 2 


d 2 x _aCL\ sin(P-0)-ato 2 cos(P-0)-Z?(O2 
dt 2 cosP 

From this equation, the linear acceleration of the slider (a s ) may be determined. 


(xix) 


The angular acceleration of the connecting rod BC (i.e. a 2 ) may be determined from 
equation (xvi) and it is given by. 


a 2 = 


a(CL\ cosB-coj 2 sin 0) — Z? CC^ sinp 
b cos P 


... {xx) 


25.5. Programme for a Slider Crank Mechanism 

The following is a programme in Fortran to find the velocity and acceleration in a slider 
crank mechanism. 


c PROGRAM TO FIND THE VELOCITY AND ACCELERATION IN A SLIDER 

c CRANK MECHANISM 

READ (*, *) A, B, E, VA, ACC, THA 
PI = 4 * ATAN (1.) 

TH = 0 

IH = 180/THA 

DTH = PI / IH 

DO 10 I = 1, 2 * I H 

TH = (I - 1) * DTH 

BET = ASIN (E - A * SIN (TH) ) / B) 

VS = - A * VA * SIN (TH - BET) / (COS (BET) * 1000) 

VB = - A * VA * COS (TH) / B * COS (BET) 

AC1 = A * ACC * SIN (BET - TH) - B * VB * * 2 
AC2 = A * VA * * 2 * COS (BET - TH) 

ACS = (AC1 - AC2) / (COS (BET) * 1000) 

AC3 = A * ACC * COS (TH) - A * VA * * 2 * SIN (TH) 

AC4 = B * VB * * 2 * SIN (BET) 

ACB = - (AC3 - AC4) / (B * COS (BET) ) 

I F (i . EQ . 1) WRITE (*, 9) 

9 FORMAT (3X,’ TH', 5X,’ BET', 4X,’ VS,’ 4X,’ VB,’ 4X,’ ACS’, 4X,’ ACB’) 

10 WRITE (*, 8) TH * 180 / P I , BET * 180 / P I, VS, VB, ACS, ACB 

8 FORMAT (6 F 8 . 2) 

STOP 

END 
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The input variables are : 

A, B, E - Length of crank AB (a), connecting rod BC (b) and offset (e) in mm, 

VA = Angular velocity of crank AB (input link) in rad/s, 

ACC = Angular acceleration of the crank AB (input link) in rad/s 2 , and 
THA = Interval of the input angle in degrees. 

The output variables are : 

THA = Angular displacement of the crank or input link AB in degrees, 

BET = Angular displacement of the connecting rod BC in degrees, 

VS = Linear velocity of the slider in m/s, 

VB = Angular velocity of the crank or input link AB in rad/s, 

ACS = Linear acceleration of the slider in m/s 2 , and 

ACB = Angular acceleration of the crank or input link AB in rad/s 2 . 

Example 25.2. In a slider crank mechanism, the crank AB = 200 mm and the connecting 
rod BC = 750 mm. The line of stroke of the slider is offset by a perpendicular distance of 50 mm. 
If the crank rotates at an angular speed of 20 rad/s and angular acceleration of 10 rad/s 2 , find at 
an interval of 30° of the crank, 1. the linear velocity and acceleration of the slider, and 2. the 
angular velocity and acceleration of the connecting rod. 

Solution. 


Given input : 


A = 200, 

OUTPUT : 

B = 750, 

E = 50, 

VA = 20, 

ACC = 10, 

THA = 30 

T H 

BET 

V S 

V B 

ACS 

AC B 

.00 

3.82 

.27 

- 5.32 

- 101.15 

- .78 

30.00 

- 3.82 

- 2.23 

- 4.61 

- 83.69 

49.72 

60.00 

- 9.46 

- 3.80 

- 2.63 

- 35.62 

91.14 

90.00 

- 11.54 

- 4.00 

.00 

14.33 

108.87 

120.00 

- 9.46 

- 3.13 

2.63 

44.71 

93.85 

150.00 

- 3.82 

- 1.77 

4.61 

55.11 

54.35 

180.00 

3.82 

- .27 

5.32 

58.58 

4.56 

210.00 

11.54 

1.29 

4.53 

62.42 

- 47.90 

240.00 

17.31 

2.84 

2.55 

57.93 

- 93.34 

270.00 

19.47 

4.00 

.00 

30.28 

- 113.14 

300.00 

17.31 

4.09 

- 2.55 

- 21.45 

- 96.14 

330.00 

11.54 

2.71 

- 4.53 

- 75.44 

- 52.61 


25.6. Coupler Curves 

It is often desired to have a mechanism to guide a point along a specified path. The path 
generated by a point on the coupler link is known as a coupler curve and the generating point is 
called a coupler point (also known as tracer point). The straight line mechanisms as discussed in 
chapter 9 (Art. 9.3) are the examples of the use of coupler curves. In this article, we shall discuss 
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the method of determining the co-ordinates of the cou¬ 
pler point in case of a four bar mechanism and a slider 
crank mechanism. 

1. Four bar mechanism 

Consider a four bar mechanism ABCD with an 
offset coupler point E on the coupler link BC, as shown 
in Fig. 25.3. Let the point E makes an angle a with BC 
in the anticlockwise direction and its co-ordinates are E 
(% >'f> 

First of all, let us find the value of BD, 7 and p . 
From right angled triangle BB { D , 

a sin 9 


and 


tan y = 


BB { 

BB X 

B X D 

AD-AB X 


[ a sin 0 

tan 1 

yd-a cos 0 


d-a cos 0 


(BD) 2 = (BB X ) 2 +(B X D) 2 = (BB X ) 2 + (AD-AB X ) 2 
= (asin 0) 2 +(d-a cos0) 2 
= a 2 sin 2 Q + d 2 + a 2 cos 2 0-2adcos0 



Fig. 25.3. Four bar mechainsm 
with a coupler point. 


= a ( sin“ 0 + cos Q) + d -2c;dcos0 

= a 2 +d 2 -2adcosQ 
Now in triangle DBC, 

,2 


cos(y + p) = 


(BD) +(BC) -(CD) 
2 BC x BD 


(cosine law of triangle) 


and 


f 2 +b 2 -c 2 


2b f 


Y+P = cos 


P = cos 


r+b 2 -c 


2 Y 


2b f 

f A +b 2 -S 


2b f 


-Y 


Let us now find the co-ordinates x E and y E . From Fig. 25.3, we find that 
x E = AE 2 = AB 1 + B x E 2 = ABj + BE X 
= a cos 9 + e cos (a + P) 

Ye = E 2 E = E 2 E\ + E X E = B X B + E X E 
= a sin 0 + e sin (a + P) 


■ ■ (i) 


(v B { E 2 = BEy ) 
... (ii) 

(v E 1 E l =B l B) 
... (iii) 


From the above equations, the co-ordinates of the point E may be determined if a, e, 9 , a 
and P are known. 
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2. Slider crank mechanism 

Consider a slider crank mechanism with an offset coupler point E, as shown in Fig. 25.4. 
Let the point E makes an angle a with BC in the anticlockwise direction and its co-ordinates are 
E (% }' F ). 

First of all, let us find the angle (3 . From right angled triangle BC X C, 

. _ BC\ BBi~BiC\ asinh-gj 


|3 =sin 


BC BC 
_l f a sin0 —q 


... (IV) 


Now 


and 



Xj : — A /'.j — \ /)] + />] /'.j — . \ />| + /> />2 

— II COS (-) I (' COS ((X |)) ... (v) 

y E = E^E = E\B 2 + B 2 E = B^B + B 2 E 

= asin9 + esin(a-(3) ... (vi) 

From the above equations, the co-ordinates of the 
point E may be determined, if a, b, e, e v 0, a and (3 
are known. 

Note : When the slider lies on the X-axis, i.e. the line of stroke of the slider passes through the axis of 
rotation of the crank, then eccentricity e t = 0. In such a case equation (iv) may be written as 


Fig. 25.4 Slider crank mechanism with 
coupler point. 


P = sin -1 


a sin 9 


25.7. Synthesis of Mechanisms 



In the previous articles, we have discussed the computer-aided analysis of mechanisms, i.e. 
the determination of displacement, velocity 
and acceleration for the given proportions 
of the mechanism. The synthesis is the 
opposite of analysis. The synthesis of 
mechanism is the design or creation of a 
mechanism to produce a desired output 
motion for a given input motion. In other 
words, the synthesis of mechanism deals 
with the determination of proportions of a 
mechanism for the given input and output 
motion. We have already discussed the 
application of synthesis in designing a cam 
(Chapter 20) to give follower a known 
motion from the displacement diagram and 
in the determination of number of teeth on 
the members in a gear train (Chapter 13) 
to produce a desired velocity ratio. 

In the application of synthesis, to 

. r , . ., ,. Note : This picture is given as additional information and is 

the design of a mechanism, the problem not a direct examp | e of the current chapter. 

divides itself into the following three parts: 



f x 
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1. Type synthesis, i.e. the type of mechanism to be used, 

2. Number synthesis, i.e. the number of links and the number of joints needed to produce 
the required motion, and 

3. Dimensional synthesis, i.e. the proportions or lengths of the links necessary to satisfy 
the required motion characteristics. 

In designing a mechanism, one factor that must be kept in mind is that of the accuracy 
required of the mechanism. Sometimes, it is possible to design a mechanism that will theoretically 
generate a given motion. The difference between the desired motion and the actual motion produced 
is known as structural error. In addition to this, there are errors due to manufacture. The error 
resulting from tolerances in the length of links and bearing clearances is known as mechanical 
error. 

25.8. Classifications of Synthesis Problem 

The problems in synthesis can be placed in one of the following three categories : 

1. Function generation ; 2. Path generation ; and 3. Body guidance. 

These are discussed as follows : 

1. Function generation. The major classification of the synthesis problems that arises in 
the design of links in a mechanism is a function generation. In designing a mechanism, the frequent 
requirement is that the output link should either rotate, oscillate or reciprocate according to a 
specified function of time or function of the motion of input link. This is known as function genera¬ 
tion. A simple example is that of designing a four bar mechanism to generate the function y = fix). 
In this case, x represents the motion of the input link and the mechanism is to be designed so that 
the motion of the output link approximates the function y. 

Note : The common mechanism used for function generation is that of a cam and a follower in which the 
angular displacement of the follower is specified as a function of the angle of rotation of the cam. The 
synthesis problem is to find the shape of the cam surface for the given follower displacements. 

2. Path generation. In a path generation, the mechanism is required to guide a point (called 
a tracer point or coupler point) along a path having a prescribed shape. The common requirements 
are that a portion of the path be a circular arc, elliptical or a straight line. 

3. Body guidance. In body guidance, both the position of a point within a moving body 
and the angular displacement of the body are specified. The problem may be a simple translation 
or a combination of translation and rotation. 


25.9. Precision Points for Function Generation 


In designing a mechanism to generate a particular function, it is usually impossible to 
accurately produce the function at more than a few points. The points at which the generated and 
desired functions agree are known as precision points or accuracy points and must be located so 
as to minimise the error generated between these points. 

The best spacing of the precision points, for the first trial, is called Chebychev spacing. 
According to Freudenstein and Sandor, the Chebychev spacing for n points in the range x s < x < x F 
(i.e. when x varies between x s and x F ) is given by 


Xj =^(x s +x F )-^(x F -x s )cos 




= — (jce +Xp)-—xAxxcos 
2 2 


2 n 
K(2j-\) 


.. (i) 



2 n 
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where x- = Precision points 


A x = Range in x = x F - x s , and 

j = 1, 2, ... n 

The subscripts s and F indicate start and finish positions respectively. 

The precision or accuracy points may be easily obtained by using the graphical method as 
discussed below. 

1. Draw a circle of diameter equal to the range A x = Xp — x§ . 

2. Inscribe a regular polygon having the number of sides equal to twice the number of 
precision points required, i.e. for three precision points, draw a regular hexagon inside the circle, 
as shown in Fig. 25.5. 

3. Draw perpendiculars from each corner which intersect the diagonal of a circle at preci¬ 
sion points x v x 2 , x 3 . 

Now for the range 1 < x < 3, x s = 1; x F = 3, and 


or radius of circle. 


A x = x F - x s = 3-1 = 2 
r = Ax/2 = 2/2 = 1 



Xo = Xc + r = Xc H-= 1-1— = 2 

2 s s 2 2 


x'i = x 2 ~ r cos 30° = x 2 



2 


and 


= 2-—cos 30° = 1.134 
2 

x 3 = x 2 + r cos 30° = x 2 + 


Fig. 25.5. Graphical method for 
cos 30° determining three precision 
points. 


Ax= x F - x s 


2 

= 2 + — cos 30° = 2.866 
2 


25.10. Angle Relationships for Function Generation 


c. 


r. 


ft 


x. 


C; 



B, 


x 


(a) Four bar mechanism. 


Fig. 25.6 


(b) Linear relationship between x and 0. 
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Consider a four bar mechanism, as shown in Fig. 25.6 (a) arranged to generate a function 
y = f(x) over a limited range. Let the range in x is (x F - x s ) and the corresponding range in 0 is 
(0 F — 0 S ). Similarly, let the range in y is ((y F - y s ) and the corresponding range in <|) is (<j)p — (j) s ). 
The linear relationship between x and 0 is shown in Fig. 25.6 (h). From the figure, we find 

that 

0 = 0s +~——(x-x s ) ••• (i) 

x F — x s 

Similarly, the linear relationship between y and (]) may be written as 


<t> = <l>s + ^—'—(y-y s ) 

yp ys 


... (ii) 



An automatic filling and sealing machine. 

Note :This picture is given as additional information and is not a direct example of the current chapter. 


For n points in the range, the equation (i) and (ii) may be written as 

0 =0 s + 9p 9s (x - -x s ) =0 s +^-(x^-x s ) 
x F — x s Ax 


and 

where 


4» / = 4>s + ^^ — As) =<\>s + / ^-(yj-ys) 

y F -y s A y 

j = 1, 2, ... n, 

A x — Xp A0 — 0p , 

Ay = - y s ’ anc ^ A <|) = (|) F - 
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Example 25.3. A four bar mechanism is to be designed, by using three precision points, to 
generate the function 

y = x 1 ' 5 , for the range l<x<4 ■ 

Assuming 30° starting position and 120° finishing position for the input link and 90° 
starting position and 180° finishing position for the output link, find the values of x, y, 9 and <|) 
corresponding to the three precision points. 

Solution : Given : x s = 1 ; x f = 4 ; 0 S = 30° ; 0 F = 120° ; <j) s = 90° ; (|)p = 180° 


Values of x 

The three values of x corresponding to three precision points (i.e. for n = 3) according to 
Chebychev’s spacing are given by 


and 


Xj =^-(x s +x F )-^(x F -x s )cos 


x l =^(1 + 4) —^(4 —l)cos 
2 2 

7t(2xl-l) 

2x3 

x 2 = ~ (1 + 4) — (4 — 1) cos 
2 2 

“ 7t(2x2-l) 

2x3 

x 3 =^(l + 4)-^(4-l)cos 

~7t(2x3-l) 

2x3 


, where j = 

1.2 Ans. 

= 2.5 Ans. 

: 3.8 Ans. 


, 2 and 3 

- (v ;= i) 

... (v j = 2) 

... (v 7 = 3) 


Note : The three precision points x v x, and x 3 may be determined graphically as discussed in Art. 25.9. 


Values of y 

Since y = x 15 , therefore the corresponding values of y are 
y l =(x 1 ) L5 = (1.2) 1 ' 5 = 1.316 Ans. 
y 2 =(x 2 ) 15 = (2.5) 15 = 3.952 Ans. 
y 3 = (x 3 ) 15 = (3.8) 1 ' 5 = 7.41 Ans. 

Note: y s =(x s ) 1 ' 5 =(1) 1-5 =1 and y F = (x F ) 15 = (4)‘ 5 = 8 

Values of 0 

The three values of 9 corresponding to three precision points are given by 


0 — 0 

= 3 ( x j ~ x s), where j - 1,2 and 3 

X F ~ X S 

120-30 

9l=30+ 4-1 {h2 ~ 1) =36 ° Ans - 

120-30 

0 2 =30+ i (2.5-1) =75° Ans. 

4-1 

120-30 

0 =30 +-(3.8-1) = 114° Ans. 

4-1 



and 
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Values of <)> 

The three values of <|) corresponding to three precision points are given by 

Jf-Js 

d). = 90+ 180-90 (1,316-1) = 94.06° Ans. 

1 8-1 

=90+ 18 ^° (3.952-1) = 127.95° Ans. 

o — l 

and (j),=90+l^——(7.41-1) = 172.41° Ans. 

3 8-1 


25.11. Graphical Synthesis of Four Bar Mechanism 

The synthesis of four bar mechanism consists of determining the dimensions of the links in 
which the output link is to occupy three specified positions corresponding to the three given positions 
of the input link. Fig. 25.7 shows the layout of a four bar mechanism in which the starting angle of 
the input link A B (link 2) of known length is 0 . Let 0 12 , 0 2 3 and 0 13 be the angles between the 
positions B 1 Z? 1 , BJ3 3 and B ] B. measured anticlockwise. Let the output link DC l (link 4) passes 
through the desired positions C v C 2 and C 3 and (J) ]2 , <t> 2 3 and <)> 13 are the corresponding angles 
between the positions C 3 C 2 , C 2 C 3 and C 3 C 3 . The length of the fixed link (link 1) is also known. 
Now we are required to determine the lengths of links B l C l and DC X (i.e. links 3 and 4) and the 
starting position of link 4 (()> ). 

The easiest way to solve the problem is based on inverting the mechanism on link 4. The 
procedure is discused as follows : 

1. Draw AD equal to the known length of fixed link, as shown in Fig. 25.8. 

2. At A, draw the input link 1 in its three specified angular positions AB V AB 2 and AB y 

3. Since we have to invert the mechanism on link 4, therefore draw a line B~D and and rotate 
it clockwise (in a direction opposite to the direction in which link 1 rotates) through an 
angle (f)| 2 (i.e. the angle of the output link 4 between the first and second position) in 
order to locate the point B' 2 . 



4. Similarly, draw another line B 3 D and rotate it clockwise through an angle (|)] 3 (i.e. angle 
of the output link between the first and third position) in order to locate point B\ . 
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B3 



Fig. 25.8. Design of four bar mechanism (Three point synthesis). 

5. Since the mechanism is to be inverted on the first design position, therefore B j and B[ are 
coincident. 

6. Draw the perpendicular bisectors of the lines B[ B\ and B' 2 B\. These bisectors intersect 
at point Cj. 

7. Join B[ Cj and C { D . The figure AB\ C X D is the required four bar mechanism. Now the 
length of the link 3 and length of the link 4 and its starting position ((f)) are determined. 

25.12. Graphical Synthesis of Slider Crank Mechanism 

Consider a slider crank mechanism for which the three positions of the crank AB {i.e. 
Bp 0 2 and 0 3 ) and corresponding three positions of the slider C {i.e. ,v |5 s 2 and ,v 3 ) are known, as 
shown in Fig. 25.9. 

In order to synthesis such a mechanism, the following procedure is adopted. 

1. First of all, draw the crank 4/J in its initial position. If the length of crank is not speci¬ 
fied, it may be assumed. 

2. Now find the ' relative poles P n and P ]3 as shown in Fig. 25.10. The relative poles are 
obtained by fixing the link A and observing the motion of the crank ,4/1 in the reverse 
direction. Thus, to find P n , draw angle YAP n equal to half of the angle between the first 
and second position (0 12 ) in the reverse direction and from AY draw IP V equal to half of 
the slider displacement between the first and second position {i.e. s 12 ). Similarly P u may 
be obtained. 

3. From P n and P l3 , draw two lines P l2 Q n and P u Q l3 such that ZAP l2 I = ZB { P l2 Qn 
and ZAP l3 1 = ZB t P n Q l3 . The lines P l2 Q l2 and P l3 Q 13 intersect at C\, which is the 
location of the slider at its first position. Now the length of the connecting rod B t C t and 
the offset (e) may be determined. 


The relative pole is the centre of rotation of the connecting rod relative to the crank rotation and 
the corresponding slider displacement. 
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(a) Three positions of the crank. 

Fig. 25.9 


3 

4-H— 

S 3 


2 


1 


s 2 


Si 2 


Si 


H - s 13 -W 


( b ) Three positions of the slider. 



Fig. 25.10 


25.13. Computer Aided (Analytical) Synthesis of Four Bar Mechanism 




(a) Four bar mechanism. ( b ) Three positions of input and output link. 


Fig. 25.11 

Consider a four bar mechanism as shown in Fig. 25.11. 

The synthesis of a four bar mechanism, when input and output angles are specified, is 
discussed below : 


Let the three positions i.e. angular displacements (0 I? 0 2 and 0 3 ) of the input link AB and 
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the three positions ((f)j, <f> 2 and <|>3 ) of the output link, as shown in Fig. 25.11 (/?), are known and 
we have to determine the dimensions a, b, c and d of the four bar mechanism. 

We have discussed in Art. 25.2 that the Freudenstein’s equation is 


k\ cos(])-/:2 cos 0 + &3 = cos(0 — (f>) 


(i) 


, d , d , , a -b +c +d z 

where h = ~ > k 2 = ’ anc ' ^3 = 


lac 


... (it) 


For the three different positions of the mechanism, the equation (t) may be written as 

k x cos 0J - k 2 cos Qj + = cos (0i - (j^) ... (iii) 

k\ cos(|)2 -k 2 cos0 2 + ^3 = cos (0 2 — § 2 ) ■■■ (iv) 

and fcj cos (|)3 — ^ cos 0 3 + ^3 = cos (03 — (|>3) ... (v) 



An off-shore oil well. 

Note : This picture is given as additional information and is not a direct example of the current chapter. 


The equations (iii), (tv) and (v) are three simultaneous equations and may be solved for k v 
k 2 and T, either by elimination method (See Examples 25.4 and 25.5) or by using Cramer’s rule of 
determinants as discussed below : 


COS^ 

COS 0! 

1 

COS(])2 

cos 0 2 

1 

COS ((>3 

cos 0 3 

1 

COS(0! 

-4>i) 

cost 

cos(0 2 

— ^ 2 ) 

cos0 

cos(0 3 

-§ 3 ) 

cosS 
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^2 — 


COS^ 
COS(f> 2 
COS (|>3 


cosCBj-^j) 1 
cos(0 2 -(|) 2 ) 1 

cos(0 3 -(|) 3 ) 1 


cos (|)| cos 0j coslBj — (f)]^) 
A 3 = cos (|>2 cos0 2 cos(0 2 — 4^2 ) 
COS(])3 COS0 3 cos(0 3 -(|) 3 ) 
Now the values of k v k , and k 3 are given by 


, A, , A 2 , A 3 
k x = —, k 2 = —— and k 3 = — 

A A A 

Once the values of k v k 2 and k } are known, then the link lengths a, b, c and d are determined 
by using equation (ii). In actual practice, either the value of a or cl is assumed to be unity to get the 
proportionate values of other links. 

Note : The designed mechanism may not satisfy the input and output angle co-ordination at positions other 
than these three positions. It is observed that a four bar mechanism can be designed precisely for five positions 
of the input and output links provided 0 and (]) are measured from some arbitrary reference rather than from 
the reference fixed link AD. In such cases, the synthesis equations become non-linear and some other means 
are required to solve such synthesis equations. 


25.14. Programme to Co-ordinate the Angular Displacement of the Input 
and Output Links 

The following is the programme in Fortran to co-ordinate the angular displacements of the 
input and output links. 

C PROGRAM TO COORDINATE ANGULAR DISPLACEMENTS OF 

C THE INPUT AND OUTPUT LINKS IN THREE POSITIONS 

READ (*, *) Ql, Q2, Q3, PI, P2, P3 
RAD = 4 * ATAN (1.0) / 180 
QA = COS (Ql * RAD) 

QB = COS (Q2 * RAD) 

QC = COS (Q3 * RAD) 

PA = COS (PI * RAD) 

PB = COS (P2 * RAD) 

PC = COS (P3 * RAD) 

AA = COS ( (Ql - PI) * RAD ) 

BB = COS ( (Q2 - P2) * RAD ) 

CC = COS ( (Q3 - P3) * RAD ) 

D = PA * (QB - QC) + QA * (PC - PB) + (PB * QC - PC * QB) 

D1 = AA * (QB - QC) + QA * (CC - BB) + (BB * QC - CC * QB) 

D2 = PA * (BB - CC) + AA * (PC - PB) + (PB * CC - PC * BB) 

D3 = PA * (QB * CC - QC * BB) + QA * (BB * PC - CC * PB) + AA * (PB * QC - PC * QB) 
A1 = D/Dl 

A2 = SQRT (A1 * A1 + A3 * A3 + 1.0 - 2 * A1 * A3 * D3 / D) 

A3 = - D/D2 

WRITE (*, 1) Al, A2, A3, 1 

1 FORMAT (6X, Al’, 7X,’ A2\ 7X,’ A3’ 7X,’ A4,’ / 4F8 . 2) 

STOP 

END 
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The input variables are : 

Q v Qt, Q 3 = Angular displacement of the 
input link AB in degrees, 

P p P^, P.. = Angular displacement of the 
output link DC in degrees. 

The output variables are : 

A, B, C, D = Ratio of length of the links AB, 
BC, CD and AD respectively 
Example 25.4. Design a four bar mechanism 
to co-ordinate the input and output angles as follows : 

Input angles = 15°, 30° and 45° ; Output angles 
= 30°, 40° and 55°. 

Solution. Given : 0j =15° ; 0 2 =30° ; 
e 3 = 45° ; (j)j = 30° ; <j) 2 = 40° ; <j) 3 = 55° 

The Freudenstein’s equation for the first position 
of the input and output link (i.e. when 0j =15° and 



(j^ = 30° ) may be written as 

k\ cos 30° - k 2 cos 15° + k 3 = cos (15° - 30°) 
or 0.866 k x - 0.966 G, + k 3 - 0.966 


Grinding machine. 

Note : This picture is given as additional 
information and is not a direct example of the 
current chapter. 

- (i) 


Similarly, for the second position (i.e. when 0 2 = 30° and (]) 2 = 40° ), 


k\ cos 40° - k 2 cos 30° + k 3 = cos (30° - 40°) 

or 0.766 k t - 0.866 k 2 + k 3 = 0.985 ... (li) 

and for the third position (i.e. when 0 3 = 45° and (]) 3 = 55° ), 


kj cos55°-k 2 cos45° + fc 3 =cos(45°-55°) 

or 0.574 k t - 0.707 k 2 + k 3 = 0.985 ... (iii) 

Solving the three simultaneous equations (/), (ii) and (iii), we get 
k x = 0.905 ; k 2 = 1.01 and k 3 =1.158 



Assuming the length of one of the links, say a as one unit, we get the length of the other 
links as follows : 

We know that k x - d/a or d = k x a = 0.905 units Ans. 














Chapter 25 : Computer Aided Analysis and Synthesis of Mechanisms • 1029 


C = d/c or c = d/C = 0.905 / 1.01 = 0.896 units Ans. 


and 


h = 


a -b +c +d z 
lac 


or -b 2 = k 3 x2ac-(a 2 +c 2 +d 2 ) 

= 1.158 x 2 x 1 x 0.896 - [l 2 + (0.896) 2 + (0.905) 2 ] 

= 2.075 - 2.622 = - 0.547 or b = 0.74 units Ans. 

The designed mechanism with AB = a = 1 unit, BC = b = 0.74 units ; CD - c - 0.896 units 
and AD = d - 0.905 units, is shown in Fig. 25.12. 

Example 25.5. Determine the proportions of four bar mechanism, by using three precision 
points, to generate y = x 1 ' 5 , where x varies between 1 and 4. Assume 9 S = 30°; A9 — 90° ; 
<f»g = 90°; and A(J> = 90° . Take length of the fixed link AD as 25 mm. 

Solution. Given : x s = 1 ; x F = 4 ; 9 S =30° ; A 9 = 9 P - 9 S = 90° ; tf> s = 90° ; 
A <t> = p F - Ps = 90° ; d — 25 mm 

We have already calculated the three values of x and y for the above given data in Example 
25.3. These values are : 

x, = 1.2 ; x 2 - 2.5 ; and x 3 = 3.8 

y, = 1.316; y 2 = 3.952 ; and y 3 = 7.41 

The corresponding values of 9 and p are 

0j = 36° ; 9 2 = 75° ; and 9 3 = 114° 
p, =94.06° ; p 2 =127.95° ; and p 3 =172.41° 

We know that the Freudenstein’s equation is 

kj cosp —k 2 cos9+ k 3 = cos(9-p) (*') 


where 


=- 


, d 

ki = — 


and 


k-x — 


—o +c + d z 
lac 


(») 


Now for the three different positions of the mechanism, the equation (i) may be written 
three times as follows : 


k { cos 94.06° - k 2 cos 36° + kj = cos (36° - 94.06°) 
or - 0.0708 k { - 0.809 k 2 + k 3 = 0.529 

Similarly k\ cos 127.95° - k 2 cos 75° + k^= cos (75° -127.95°) 
or - 0.615 k x - 0.259 k 2 + k 3 = 0.6025 

and k x cos 172.41° - k 2 cosl 14° + k 2 = cos (114°-172.41°) 

- 0.9912 k x + 0.4067 k 2 + k 3 = 0.5238 
Solving three simultaneous equations (iii), ( iv ) and (v), we get 
k x = 0.6 ; k 2 = 0.453 ; and k 3 = 0.12 
Now from equation (ii), 


.. (iii) 
■ ■ (iv) 

... (v) 



— = 41.7mm 
0.6 


Ans. 
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and 


d 25 . 

c = — =-= 55.2nim Ans. 

k 2 0.453 


< = (a^+c 2 +d — k 2 Y.2ac)^ 


= [(41.7) 2 + (55.2) 2 +(25) 2 -0 


-0.12x2x41.7x55.2 


: 69.7 mm Ans. 


The designed four bar mechanism AB^ C 2 D in one position ( i.e. for 0 2 , x 2 and (j) 2 , y 2 ) is 
shown by thick lines in Fig. 25.13. 

10 1.316 



The other two positions of the four bar mechanism may be drawn by joining I) , C 2 (i.e. 
0j, x 1 and <])], y;) and B 3 C 3 (i.e. © 3 , *3 and < 1 ) 3 , y 3 ). 

Note : In the above example, the motion of input link and output link is taken clockwise. 

Example 25.6. Synthesize a four bar linkage, as shown in Fig. 25.14, using Freudenstein's 
equation to satisfy in one of its positions. The specification 
of position 6, velocity CO and acceleration a are as follows : 

0 = 60°, © 2 = 5 rad/s; a, = 2 rad/s 2 ; 

(b = 90°; ©, = 2 rad/s; a, = 7 rad/s 2 . 


Solution : Given : 0 = 60° ; ©., = 5 rad/s ; a, = 2 rad/s"; 
= 90°, © 4 = 2 rad/s; a 4 = 7 rad/s 2 

The four bar linkages is shown in Fig. 25.15. Let 

AB = Input link = a, 

BC = Coupler = b, 

CD - Output link = c, and 

AD = Fixed link = d. 

The Freudenstein’s equation is given by 

costj)- k 2 cos 0 + k 2 =cos(0-<j>) ... (i) 




Fig. 25.15 
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where 


k\ - —; k 2 — —; and *3 = 


a-b z +c z + d z 
2 ac 


Substituting the value of 0 and <f> in equation (i), 

*2 cos 90° — * 2 cos 60° + * 3 = cos(60° — 90°) 

*2 x 0-* 2 xO .5 + ^3 = 0.866 

-0.5 * 2 + * 3 =0.866 

Differentiating equation (/) with respect to time, 

. d§ . d 0 . d (0 —(j)) 

*2 x-sm<|)x — -k 2 x-sin 0 x— = -sin( 0 -<j))x-— 

dt dt dt 

-k\ sin (f)t 0 4 +k 2 sin 0 ( 0 2 = -sin (0 - tJ>Xco 2 - (O 4 ) 


(») 


{Hi) 


—*2 x sin 90° x 2 + k 2 sin 60° x 5 = — sin(60°—90°) (5 — 2) 

5V3 _ 3 

— 2k] H- kn — — 

2 2 

*2 =2.165 k 2 -0.75 

Now differentiating equation (iii) with respect to time. 


~h 


, d( 0 4 c/(|) 

sm(])X—- + co 4 cos(])X— 1 
dt dt 


+ kn 


<r/(b , r/0 

=C 0 4 ; and—= 
dt dt 


■ ■ (iv) 


. dec, c/0 

sin0x—- + C 0 i cos0x — 
dt ‘ dt 


or 


• /a 1, d(® 2 — CO4) , .. ,, c/(0-<))) 

= -sin(0-()))---—+ (co 2 -C 0 4 ) cos(0-<j))x- 


dt 


dt 


-k\ rsin<|)xa 4 +CO 4 cos(f)J + * 2 |^sin 0 xa 2 +(Oj cos ©J 

= — [^sin( 0 —(|>) (a 2 -a 4 )+(G ) 2 — co 4 ) 2 cos( 0 — 4 >)j| 

-k x j^sin 90 ° x 7 + 2 2 cos 90 °J+ k 2 j^sin 60 ° x 2 + 5 2 cos 60 °J 

= ~[sin( 60 ° - 90 °) ( 2 — 7 ) + (5 — 2 ) 2 cos( 60 ° - 90 °)] 

-k x (7 + 0) + k 2 (1.732 +12.5) = - (2.5 + 7.794) 

-7* t +14.232* 2 =-10.294 

* t = 2.033 * 2 + 1.47 ... (v) 

From equations (tv) and (v), 

2.165*2-0.75 = 2.033 * 2 + 1.47 or * 2 = 16.8 













1032 • Theory of Machines 


From equation (v) 

k x = 2.165&2 -0.75 = 2.165x16.8-0.75 = 35.6 

and from equation (ii), 

k 3 = 0.5 k 2 + 0.866 = 0.5 x 16.8 + 0.866 = 9.266 
Assuming the length of one of links say a as one unit, we get the length of the links as 
follows : 

We know that k^ = d / a or d = k x . a - 35.6 units Ans. 

k 2 = d / c or c - d / k 2 = 35.6 / 16.8 = 2.12 units Ans. 


and 


a 2 -b 2 +c 2 +d 2 _ l 2 -b 2 +(2.12) 2 +(35.6) 2 
2ac 2x1x2.12 



l-b 2 +4.494 + 1267.36 1272.854-6 2 

y. zoo =-=- 

4.24 4.24 

b 2 = 1272.854 - 9.266 x 4.24 = 1233.566 
b — 35.12 units Ans. 

Example 25.7. Synthesize a four-bar mechanism to generate a function y = sin x for 
0 < x < 90°. The range of the output crank may be chosen as 60° while that of inut crank be 120°. 
Assume three precision points which are to be obtained from Chebyshev spacing. Assume fixed 
link to be 52.5 mm long and 9j =105° and (j)| = 66°. 

Solution. Given : ^s=0; x p =90°; A(f) = 60°; A0 = 12O°; ii = 52.5mm; 

0! =1.05°; (j) 1= 66° 

The three values of x corresponding to three precision points (i.e. for n = 3), according to 
Chebyshev spacing are given by 


and 


Xj = Os +-%)--U F -x s )cos 
2 

x, =-(0 + 90)--(90-0) cos 
2 2 

= 45 -45 cos 30° = 6° 

x 2 =^(0 + 90)-^(90-0) cos 

= 45-45 cos 90° = 45° 

x =i(0 + 90)--(90-0)cos 
2 2 

= 45-45 cos 150° = 84° 

Since y = sin x, therefore corresponding values of y are 

yi = sin jq = sin 6 ° = 0.1045 


71(2x3-1) 

2x3 


71(2x2-1) 

2x3 


n(2j-l) 
2 n 

71(2x1-1) 

2x3 


where j = 1, 2, 3 


...(v 7=1) 


■■■(••■ 7 = 2) 
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and 


and 


and 


and 


y 2 = sin x 2 = sin 45° = 0.707 
y 3 = sin x 3 = sin 84° = 0.9945 
y s = sin x s = sin 0° = 0 
y F = sin x F = sin 90° = 1 

The relation between the input angle (0) and x is given by 


Also 


0 , — 0c + 


0c -0c 


(xj x s ), where j = 1, 2 and 3. 


x F x s 

The above expression may be written as 
A0 

Q J =B s + —(Xj-x s ) 

Ax 

The three values of 0 corresponding to three precision points are given by 

.. (v x s = 0 )... 0 ) 
A0 


“ ®s + . Xx i 
Ax 


09 — 0 s~t x x 2 
Ax 


A0 


0 3 = 0 S -l-xx. 

Ax 

From equations (/), 07) and 077), 


02 — 01 = — (*2 - ^ 1 ) = 7 T 7 (45 - 6 ) = 52° 
Ax 90 


e 3 -0 2 = ^(x 3 -x 2 ) = ^(84-45) = 52° 
Ar 90 


03 -0! = fu 3 -i) = ~~~ (84 — 6 ) = 104° 
Ax ^ 90 


... 07) 

.... 077) 

... 0 v) 

(y Ax = Xp — x s = 90 — 0 = 90) 

... (v) 

... 0V) 


Since 0j = 105° (Given), therefore 

02=0!+ 52° = 105° + 52° = 157° 

0 3 = 0 2 +52° = 157°+ 52 = 209° 

The relation between the output angle (cf)) and y is given by 


4>y=<t> S +^-^ 
y f y s 

This expression may be written as 


(yj Ts)> when j = 1, 2 and 3 


, A<|) 

= ^>S +—(>'; -Js) 
Ay 
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The three values of (]) corresponding to three precision points are given by 


and 


4>i =^>s + —x^i 
Ay 

, , A4» 

^2 — + . X y2 

Ay 

Aty 

<t>3 =^>S + — X J3 

Ay 

From equations (vii), (viii) and ( ix ), 


, — (> x = — (y 2 — = -(0.707-0.1045) = 36.15° 

Ay 1 


.... (v y s = 0 ) ... (vii) 
... (viii) 

... (ix) 

... (x) 


(> 3 - <|> 2 = — (y 3 - y 2 ) = — (0.9945 - 0.707) = 17.25 
Ay 1 


('•• Ay = y F -y s = 1 - 0= 1) 
... (xi) 


<>3 — <,j = — (y 3 -yj) = -(0.9945-0.1045) =53.4° 
Ay 1 


.. (xii) 


Since (jjj = 66 ° (Given), therefore 

<> 2 =<)>! +36.15° = 66 ° + 36.15° = 102.15° 

(j>3 = <j ) 2 +17.25° = 102.15° + 17.25° = 119.40° 

We have calculated above the three positions i.e. the angular displacements (0 p 0 2 and 0 3 ) 
of the input crank and the three positions (t) p <J ) 0 and q,) of the output crank. Now let us find the 
dimensions of the four bar mechanism. 

Let a = Length of the input crank, 

b = Length of the coupler, 
c - Length of the output crank, and 
d = Length of the fixed crank = 52.5 mm 

We know that the Freudenstein displacement equation is 
kj cos q - k 2 cos 0 + k 3 = cos (0 - q) 


(Given) 

... (xiii) 


where 


, d, 

k \ =-; 


k-, = — and 


kn = 


a -b +c +d z 
2ac 


The equation (xiii) for the first position of input and output crank (i.e. when 0 = 45° and 
66 °) may be written as 

kj cos — k 2 cos 0 [ + k 3 = cos( 0 j - (])[) 

kj cos 66 ° — k 2 cos 105° + k 3 = cos(105° — 66 °) 

0.4067k! +0.2588k 2 +k 3 = 0.7771 


... (xiv) 
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Similarly, for the second position (i.e. when = 157° and <S? n = 102.15°), 
cost])? — k 2 cos0 2 + k 3 = cos(0 2 ~§ 2 ) 


k x cos 102.15° - k 2 cos 157° + k 3 = cos(l 57° -102.15°) 

-0.2105^ + 0.9205& 2 + k 3 =0.5757 ...(xv) 

and for the third position (i.e. when 0 3 = 209° and p, = 119.4°), 
kj cos <(>3 — k 2 cos 0 3 + k 3 = cos(0 3 - (f) 3 ) 


fcj cos 119.4° - k 2 cos 209° + k 3 = cos(209° -119.4°) 


- 0.4909&J +0.8746 k 2 +k 3 = 0.007 ... ( xvi) 

Solving the three simultaneous equations (xiv), (xv) and (xvi), we get 
k l = 1.8; k 2 = 1.375 and k 3 = - 0.311 

Since the length of the fixed link (i.e. d = 52.5 mm) is known, therefore we get the length 
of other links as follows: 

We know that 

k x - d / a or a = d / k l = 52.5 / 1.8 = 29.17 mm Ans. 

k 2 -d/c or c = dlk 2 = 52.5 / 1.375 = 38.18 mm Ans. 


and 


ki = 


a -b +c +d 
2 ac 


or b 2 - a 2 +c 2 +d 2 -k 3 x2ac 

= (29.17) 2 + (38.18) 2 + (52.5) 2 - (-0.311)x 2 x 29.17 x 38.18 = 5758 
b = 75.88 mm Ans. 



25.15. Least Square Technique 

Most of the mechanisms are not possible to design even for five positions of the input and 
output links. However, it is possible to design a mechanism to give least deviation from the specified 
positions. This is done by using least square technique as discussed below : 

We have already discussed that the Freudenstein’s equation is 
kj cos (f) - k 2 cos 0 + k 3 - cos (0 — (f)) = 0 

The angles 0 and (|) are specified for a position. If 0, and <J> ( are the angles for z'th 
position, then Freudenstein’s equation may be written as 

k, cos <]) ( - - k 2 cos 0 ; - + & 3 - cos (0 ; - - (j) ( -) = 0 

Let e be the error which is defined as 


e = cos(j) ; — k 2 cos 0 ; - + fc 3 — cos (0 Z - -()), )] 2 
i=l 

For e to be minimum, the partial derivatives of e with respect to k v k 2 , k 3 separately must 
be equal to zero, i.e. 


de_ 

dk l 


= 0 


= 0, and 
dk 2 


de_ 

dk 3 


= 0 
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= 2 V' [k\ cos §j— k 2 cos 9 ( - + k 2 - cos (0 ( - - <j) ; -)] cos (j> ; = 0 
dk, “ 


fcj^COS 2 (|) ; - -^2^COS0 ; -COS(|) ; -+^3^COS(|) ; - = ^cos(0 ( - - (|) ; -) COS (j) ; - ■■■ (0 

=1 i=1 


=1 i=1 


Similarly, jf- = -2^ [ftj cos (|) ; - - k 2 cos 0 ; + k 2 - cos (0 ; - - (|) ; -)] cos 0 ; =0 

a/: 2 f-r 


^ cos ()),■ cos 0 ; + k 2 ^ cos 2 0 ; - + k 2 ^ cos 0,- = ^ cos (0 ( - - ()),■) cos Qj ... (ii) 


3^ ' 

j^° w = 2^ [k\ cos ()),■ - cos 0 ; - + k 2 - cos (0 ; - - (j) ( -)] = 0 


"3 «=i 


or /cj^cos^y + k 2 ^^cosQ j + & 3 ^1 = ^cos(0 ( --(]),■) ... (Hi) 

i =1 i=l (=1 i=l 

The equations (i), (ii) and (iii) are simultaneous, linear, non-homogeneous equations in 
three unknowns k v k 2 and k y These equations can be solved by using Cramer’s rule. 

25.16. Programme Using Least Square Technique 

The following is a programme in Fortrans to find the ratio of lengths for different links by 
using the least square technique. 

The input variables are : 

J = Number of specified positions 

TH ( I) = Angular displacements of the input link AB for / = 1 to J (degrees), and 
PH ( I) = Angular displacements of the output link DC for I - 1 to J (degrees). 
The output variables are : 

A, B, C, D = Ratio of the lengths of the links AB, BC, CD and AD respectively. 

C PROGRAM TO COORDINATE ANGULAR DISPLACEMENT OF THE 

C INPUT AND OUTPUT LINKS IN MORE THAN THREE POSITIONS TO 

C FIND RATIO OF DIFFERENT LINKS USING LEAST SQUARE TECHNIQUE 

DIMENSION 
READ (*, *) I 

READ (*, *) (TH (I), I = 1, J), PH (I), I = 1, J) 

RAD = 4 * ATAN (1.0) / 180 
DO 10 K = 1 . I 

A1 = A1 + (COS (PH (K) * RAD ) ) * * 2 

A2 = A2 + (COS (TH (K) RAD ) ) * (COS (PH (K) * RAD ) ) 

A3 = A3 + (COS (PH (K) * RAD ) ) 

B1 = A2 

B2 = B2 + (COS (TH (K) * RAD ) ) * * 2 
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B3 = B3 + (COS (TH (K) * RAD ) ) 

PI = A3 
P2 = B3 
P3 = J 

TT = COS ( ( TH (K) - PH (K) * RAD ) 

Q1 = Q1 + TT * COS (PH (K) * RAD ) 

Q2 = Q2 + TT * COS (TH (K) * RAD 
10 Q3 = Q3 + TT 

D = At * (B2 * P3 - B3 * P2) + B1 * (P2 * A3 - P3 * A2) + PI * (A2 * B3 - A3 * B2) 

D1 = Q1 * (B2 * P3 - B3 * P2) + B1 * (P2 * Q3 - P3 * Q2) + PI * (Q2 * B3 - Q3 * B2) 

D2 = A1 * (Q2 * P3 - Q3 * P2) + Q1 * (P2 * A3 - P3 * A2) + PI * (A2 * Q3 - A3 * Q2) 

D3 = A1 * (B2 * Q3 - B3 * Q2) + B1 * (Q2 * A3 - Q3 * A2) + Q1 * (A2 * B3 - A3 * B2) 

Q = D / D1 
R = - D / D2 

P = SQRT (Q * Q + r * r + 1. - 2. * r * r * 03 / D) 

WRITE (* , 9) Q, P, r, 1. 

9 FORMAT (6X, ’ Q\ 7X,’ P\ 7X,’ r\ 7X,’ D’ / 4F8 . 2) 

STOP 

END 

25.17. Computer Aided Synthesis of Four Bar Mechanism with Coupler 
Point 

Consider a four bar mechanism ABCD with a couple point E, as shown in Fig. 25.16, 
which is specified by r and Y . 



Fig. 25.16. Four bar mechanism with a couple point. 

Let 8 ,, 0 2 and 0 3 = Three positions of the input link AB, 

r v r , and r , = Three positions of the coupler point E from point O, and 
Yl, y 2 and Y 3 = Three angular positions of the coupler point E from OX. 

The dimensions a, c, e, / and the location of points A and D specified by (q, (3 ) and (p, a ) 
respectively, may be determined as discussed below : 
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Considering the loop OABE, the horizontal and vertical components of vectors q, a, e and 


r are 

^cos(3 + acos9 + ecos8 = rcosy ...(/) 

and <2 , sin[3 + asin0 + esinS= rsiny —(it) 

Squaring equations (i) and (ii) and adding in order to eliminate angle 8, we have 

g[2rcos(y-P)] + a[2rcos(0-y)] + e 2 -q 2 -a 2 = r 2 + ya[2cos(0-|3)l ... (iii) 

Let q = itj ; a = k 2 ; e 2 -q 2 -a 2 =k 3 ; and q a - k 4 = k 1 k 2 ... (tv) 

Now the equation (iii) may be written as 

[2rcos(y-p)] + it 2 [2rcos(0 — y)] + jt 3 = r 2 +k 4 [2cos(0-p)] ... (v) 


Since k 4 = k\k 2 , therefore the equation (v) is difficult to solve for k\,k 2 ,k 3 and k 4 . Such 
type of non-linear equations can be solved easily by making them linear by some substitutions as 
given below : 

Let kj=/j+Xt;ij ; k 2 =l 2 +Xm 2 ; and k 3 =l 3 +X m 3 ... (vi) 

where A = k 4 = Zy k 2 = (/j + A, t/ij) (l 2 + A m 2 ) 

= l\ l 2 +/[ Am 2 + A«ij l 2 + A 2 tM| m 2 
or m 2 A" + (Zj m 2 + 1 2 ni\ — 1) A + l\ l 2 =0 

or AX 2 + BX + C = 0 

-B±Jb 2 -4AC 

A =--- ... (vtt) 

2A 

where A = m 3 m 2 ; B = (Zj m 2 + 1 2 m 3 -1) ; and C = l 2 ... ( viii ) 

Substituting the values of k v k 2 , Zy and Zy in equation (v), 

(/| + Atn|)[2rcos (y—P)] + (l 2 + Am 2 )[2 rcos(0 — y)] + (Z 3 + A m 3 ) 

= r 2 +A([2cos(0-p)] 

Equating the terms with A and without A separately equal to zero, we get the components 
into two groups, one with A and the other without A ■ These components are 


Z 1 [2rcos(y-p)] + Z 2 [2rcos(0-y)] + Z 3 = r 2 ... (ix) 

and nq[2rcos(y-P)] + «i2[2rcos(0 —y)] + m 3 = 2cos(0-P) ... (jc) 

The equation (ix) for the three positions of 0 , r and Y may be written three times as 
follows : 

Z 1 [2r 1 cos(y 1 -P)] + Z 2 [2r 1 cos(0 1 -y 1 )] + Z 3 =rf ... (xi) 

Z 1 [2r 2 cos(y 2 -P)] + Z 2 [2r 2 cos(0 2 -y 2 )] + Z 3 = r 2 ... (xii) 

Z 1 [2r 3 cos(y 3 -P)]+Z 2 [2r 3 cos(0 3 -y 3 )] + Z 3 = r 3 2 ... (xiii) 


Similarly, equation (jc) for the three positions of 0 , r and Y may be written three times as 
follows : 
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«i 1 [2rcos(Yi-(3)l + m 2 [2rcos(0i — Yi)] + ?n 3 =2cos(0j-(3) ... ( xiv ) 

«! 1 [2rcos(y 2 -P)] + «i 2 [2rcos(0 2 -Y 2 )] + »*3 = 2cos(0 2 ~P) ... (jcv) 

;» 1 [2rcos(Y 3 -P)] + /n 2 [2rcos(0 3 -Y 3 )] + «?3 = 2cos(0 3 ~P) ... (xvi) 

The equations ( xi ), ( xii ) and ( xiii ) are three linear equations in l v l 2 , Z 3 . Similarly, equations 
(xiv), (xv) and (xvi) are three linear equations in m v m 2 and Assuming a suitable value of p , 
the values of l v l 2 , / 3 and m v m 2 , m 3 may be determined by using elimination method or Cramer’s 
rule. 

Knowing the values of l v l 2 , Z 3 and m v m 2 , m y we can find the value of X from equation 
(vii). Now the values of k v k 2 and Zc, are determined from equation (vi) and hence q, a and e are 
known from equation (iv). Using equation (i) or (ii), we can find the three valves of 5 Z.e. 8], S 2 
and S 3 . From equation (/), we 1 
e cos 5 = r cos 

8j = cos -1 
Similarly, § 2 = cos -1 


ave 

y-^cosp-acos 0 
q cos Yi —q cos P - a cos 0 ] 


r 2 cos Y 2 —q cos P - a cos 0 2 


... (xvii) 

... (XVIII) 


and 


5 3 = cos 


r 3 cos Y 3 - q cos P - a cos 0 3 
e 


... (xix) 


Thus by considering the loop OABE, we can find the values of q, a, e, P and 5 ■ 

Now considering the loop ODCE in order to find p, c, f a and V ■ The horizontal and 
vertical components of vectors p, c, f and r are 

p cos a + c cos <j) + / cos Y|/ = rcosy ... (xx) 


and p sin a + c sin (f> + / sin \|/= r sin y ... (xxi) 

Since these equations are similar to equations (i) and (ii), therefore we shall proceed in the 
similar way as discussed for loop OABE. 

Squaring equations (xx) and (xxi) and adding in order to eliminate angle (f), we have 


p[2rcos(Y-a)] + /[2rcos(\|/-Y)] + c 2 - p 2 - f 2 = r 2 + p/[2cos(\|/-a)] ■■■ (xxii) 

Let p = k 5 ;/= k 6 ■, c 2 - p 1 -f ~ = k 2 and p f = k & = k 5 k 6 ... (xxiii) 

Now equations (xxii) may be written as 

£ 5 [2rcos(Y-a)] + A: 6 [2rcos(\i/-y)] + £ 7 = r 2 + £ 8 [2cos(\|r-a)] ••• (xxiv) 

The equation (xxiv) is a non-linear equation and can be solved easily by making it linear 
by some substitutions as given below : 

Let k 5 = Z 5 +A,j m 5 ; k 6 =l 6 +X l m 6 ; and k 2 =l-j +X l m-j ... (xxv) 

where X l = k & = k 5 k 6 = (Z 5 + A,j m 5 ) (Z 6 + A,j m 6 ) 

= h l 6+ 1 5 \ m 6 + h m 5 h + m 5 m 6 


m 5 m 6 X] + (Z 5 m 6 + Z 6 m 5 - 1) X t +1 5 1 6 = 0 



or 
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or 


DX\ + E’k l + F = Q 




-Et^E 2 -4DF 
2D 


... (xxvi) 


where D = m 5 m 6 ;E = (l 5 m 6 +1 6 m 5 —1) ; and F - l 5 l 6 ... (xxvii) 

Substituting the values of k 5 , k (] , k 7 and k K in equation (xxiv), 

(l 5 +A,j ra 5 )[2rcos(y-a)] + (/ 6 +A,j w 6 )[2rcos(\|/-y)] + / 7 +X l m 7 

= r 2 +A, 1 [2cos(\)/-a)] 

Equating the terms with A, and without A, separately equal to zero, we get the components 
into two groups, one with A, and the other without A, ■ These components are 

l 5 [2r cos (y- a)] + Z 6 [2r cos (\|/ - y)] + / 7 = r 2 ■■■ (xxviii) 

and m 5 [2rcos(y-a)] + /n 6 [2rcos(t(/-y)] + /?j 7 =2cos(\|/-a) ... (x:ux) 

The equation ( xxviii ) for the three positions of r, y and t|/ may be written three times as 
follows : 


/ 5 [2q cos (y t - a)] + / 6 [2r x cos (Vt - Yi)] + 1 7 = r 2 ... {xxx) 

Z 5 [2r 2 cos(y 2 —CX)J + /^[2/ 2 cos(v|/ 2 —y 2 )J + / 7 — r 2 ... (xxxi) 

Z 5 [2r 3 cos(y 3 —(X)J + /^[2r 3 cos(v|/ 3 —y 3 )J + / 2 = r 3 ... (xxxii) 

Similarly, equation (xxix) for the three positions of r, y and V may be written three times 
as follows : 

m 5 [2qcos(yj -a)] + m 6 [2qcosCti/j-yj)] + m 7 = 2008 ( 1 )/! -a) ... (xxxiii) 

n/ 5 [2r 2 cos(y 2 -a)] + «; 6 [2r 2 cos(\(/ 2 -y 2 )] + ;;j 7 = 2cos(\(/ 2 -a) .. (xxxiv) 

«! 5 [2r 3 cos(y 3 -a)] + «i 6 [2r 3 cos(\(/ 3 -y 3 )] + n; 7 = 2cos(\(/ 3 -a) ... (xxxv) 


The equations (xxx). (xxxi) and (xxxii) are three linear equations in / 5 , / 6 and l T Similarly, 
equations (xxxiii), (xxxiv) and (xxxv) are linear equations in m y m 6 and m 7 . Assuming a suitable 
value of a, the values of l 5 , l 6 , l 7 and m y m 6 , m 7 may be determined by using elimination method 
or Cramer’s rule. 


Knowing the values of / 5 , / 6 , l 7 and m 5 , m (y m 7 , we can find the value of A,j from equation 
(xxvi). Now the values of k y k 6 and k 7 are determined from equation (xxv) and hence p, f and c are 
known from equation (xxiii). 

Assuming the value of \|/j, the corresponding values of \j/ 2 and t|/ 3 may be calculated as 
follows : 

Since the angular displacements of the coupler link BCE is same at the points B and C, 
therefore 


¥ 2 ~¥i — S 2 — S, 

or \|/ 2 = xj/j + (8 2 - Sj) ... (xxxvi) 

Similarly, \|/ 3 = \|/j + (5 3 -8]) ... (xxxvii) 

If the mechanism is to be designed for more than three positions of the input link AB and 
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the same number of positions of the couple point E, then the least square technique is used. The 
error function from equations ( ix ) and ( x ) are defined as 

<?1 = cos(y— p)} + / 2 {2rcos (0 — y)} + Z 3 -r 2 ] 2 ... (xxxviii) 

and e 2 =^[Ht 1 {2rcos(y-p)} + m 2 {2rcos(0-y)} + »t 3 -2cos(0-p)] 2 ... (xxxix) 



An aircraft assembling plant. 

Note :This picture is given as additional information and is not a direct example of the current chapter. 


For e x and e 2 to be minimum, the partial derivatives of e x with respect to l v l 2 , Z 3 and 
partial derivatives of e 2 with respect to m x , m 2 , /;?, separately must be equal to zero, i.e. 


^L =0 ;^L=0; ^1=0 
3 /[ 3/ 2 3 /^ 


and 


^ 2 _ = 0 . de 2 


dim 


3 Win 


= 0 ; 


‘3 

de 2 

dm-. 


= 0 


First consider when 


—^ = 0 

3/. 


... (xxxx) 


or 


^2^Z 1 2rcos(y-(3) + / 2 2rcos(0-y) + /3 -r 2 J 2rcos(y-(3) = 0 
1 

[2r cos (y - p)] 2 + / 2 [2r cos (0 - y)] [2r cos (y - (3)] 

i t 

+Z 3 [2r cos (y - p)] = cos (Y'“ P)1 r 2 

1 1 




! ^ [2r cos (y - P)] [2r cos (0 - y)] + Z 2 ^ [2r cos (0 - y)] 2 
l l 

+Z 3 ^2[2rcos(0-y)] = ^[2rcos(0-y)]t 


... (xxxxi) 


... (xxxxii) 
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and 


for !r = 0 ^E[2rcos(Y-p)] + Z 2 y[2rcos(0-Y)] + / 3 yi = y/- 2 

a/ 3 T 1 T 

The above three equations can be solved by using Cramer’s rule to find l v /, and l v 


(xxxxiii) 


In the similar way as discussed above, for 


= 0 


m x E [2r cos (y -P)] 2 
1 


+ m 2 E [2r cos (0 - y)] [2r cos (y — (3)] 

1 

+m 3 Ecos(y-p)] = E[2rcos(0 — (3)][2rcos(y — (3)1 ... (xxxxiv) 
1 1 


Similarly, for 


de 2 

dm 2 


= 0, 


m j E [2r cos (y - p)] + [2r cos (0 - y)] + m 2 E [2r cos (0 — y)] 2 

1 

+m 3 E [2r cos (0 - y)] = E pr cos (0 - p)] [2r cos (0 — y)] ( yyyy v ) 

1 1 


n?jE [2rcos(y- P)] + m 2 E [2rcos (0 -y)] + m 3 E1 = E cos (0 — y)] ... ( xxxxvi ) 
1 1 1 
The above three equations can be solved by using Cramer’s rule to find m v m 2 and m y 
Knowing the values of l v l 2 , Z 3 and m v m 2 , m y we can find the value of X from equation 
(vii) and k v k 2 , k 3 from equation (vi). Thus q, a and e are determined. Now 5|,5 2 ,§ 3 may be 
determined by using equation (i) or (if). 

The values of p, c and / are obtained by solving equation (xxiv) in the similar way as 
discussed earlier. 

25.18. Synthesis of Four Bar Mechanism For Body Guidance 

Consider the three positions of a rigid planer body containing the points A and B as shown 
in Fig. 25.17 (a). The four bar mechanism for body guidance, considering the three positions of the 
body, may be designed graphically as discussed below. 

1. Consider the three positions of the points A and B such as A p A 2 , A 3 and B v B 2 , B 3 as 
shown in Fig. 25.17 (a). 

2. Find the centre of a circle which passes through three points A 3 , A 2 , A 3 . This is obtained 
by drawing the perpendicular bisectors of the line segments A l A 2 and A 2 A 3 . Let these bisectors 
intersect at 0 A . It is evident that a rigid link A 0 A pinned to the body at point A and pinned to the 


and for 




= 0, 
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ground at point 0 A will guide point A through its three positions A,, A 2 and A y 
Y 




(a) 

Fig. 25.17. Four bar mechanism for body guidance. 

3. Similarly, find the centre 0 B of a circle which passes through three points B v /A, B y It 
is evident that a rigid link B 0 B pinned to the body at point B and pinned to the ground at point 0 H 
will guide point B through its three positions B v B 2 and B y 

4. The above construction forms the four bar mechanism 0 A ABO b which guides the body 
through three specified positions. Fig. 25.17 ( b ) shows a four bar mechanism in these three positions. 
The points O A and 0 B may be determined analytically as discussed below : 

Consider the three positions of the point A such as A p A 2 , A y Let the co-ordinates of these 
points are A l (x v y; A n (x 2 , y 2 ) and A 3 (x 3 , y 3 ). Let the co-ordinates of the point 0 A are (x, y). 
Now we know that 

Distance between points Aj and 0 A , 


A 1 C A =[(x 1 -x) 2 + (y 1 -y) 2 ] 1/2 


Similarly, distance between points A n and 0 A , 
A 2 O a — [(^2 

and distance between points A 3 and O A , 


A 2 0 A =[(x 2 -x) 2 +(y 2 -y) 2 f /2 


A 3 0 A =[(x 3 -x) 2 +(y 3 -y) 2 ] 1/2 


■ ■ (i) 


(«) 


(in) 


For the point 0 A to be the centre of a circle passing through the points A p A, and A 3 , the 
distances A 3 O a , A 2 O a and A 3 O a must be equal. In other words, 

AjO a = A 2 O a = A^Oa 


Now considering AjO a = A 2 O a , we have 

-, 1/2 


[(^1 -x) 2 +(y 1 -y) 2 J =[(x 2 -M 2 + (A 2 — A) 9 ] 

Similarly, considering A 2 O a = A 3 O a , we have 

^(x 2 -x) 2 + (y 2 -y) 2 J =[(x 3 -x) 2 + (y 3 -y) 2 ] 


■ ■ (iv) 


... (v) 


Squaring both sides of the equations (tv) and (v) and simplifying, we get the following two 
equations in the unknowns x and y. 

2x (x 2 -X]) + 2y (y 2 - y,) + (x 2 -x\) + (y 2 - y 2 ) = 0 ... (vt) 

2x(x 3 — x 2 ) + 2y (y 3 — y 2 ) + (x 2 — x 3 ) + (y 2 — y 2 ) = 0 (vii) 




and 


... (vt) 
... (vii) 
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The equations (vi) and (vii) are simultaneous equations and may be solved to find the 
co-ordinates x, y of the point O A . This point 0 A becomes the location of the fixed pivot guiding the 
point A. The length of the guiding link O a A may be determined by any of the equations (i), (ii) or 
(hi). 

In the similar way, as discussed above, we can find the location of the fixed pivot point 0 B 
and the length of the link O b B. 

Example 25.8. Synthesize a four bar mechanism to guide a rod AB through three consecu¬ 
tive positions AjBj, Aff and A { /I { as shown in Fig. 25.18. 



Fig. 25.18 

Solution : In order to synthesize a four bar mechanism, we shall use the graphical method 
as discussed below : 

1. Join points A p A, and A,, A v Draw the perpendicular bisectors of line segments A X A 2 
and Ay4 3 to intersect at 0 A , as shown in Fig. 25.19. It is evident that a rigid link O a A j pinned to the 
body at point A x and pinned to the ground at point 0 A will guide point A x through its three posi¬ 
tions. 


* 
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2. Similarly, join points B v II, and B 2 , B y Draw the perpendicular bisectors of line seg¬ 
ments B JC and B 1 B 3 to intersect at 0 B as shown in Fig. 25.19. It is evident that a rigid link O r B 1 
pinned to the body at point />' and pinned to the ground at point 0 B will guide point B, through its 
three positions. 

3. From above we see that the points 0 A and 0 B are the required fixed points and 
0 A Aj B j 0 H is one position of the four bar mechanism. The other two positions of the mechanism 
will be 0 A A 2 B 2 Oq and 0 A A 3 B 2 Oq . 


25.19. Analytical Synthesis for Slider Crank Mechanism 

A slider crank mechanism is shown in Fig. 25.20. In the sysnthesis problem of the slider 
crank mechanism, the displacement (s) of the slider C has to co-ordinate with the crank angle (0 ) 
in a specified manner. For example, consider that the displacement of the slider is proportional to 
crank angle over a given interval, i.e. 


s-s 5 = C(0-0 S ), for 0 S < 0 < 0 F 


... (i) 


where C = Constant of proportionality, and 

s = Displacement of the slider when crank angle is 0 . 


The subscripts s and F denote 
starting and finishing positions. 

The synthesis of a slider crank 
mechanism for three precision points 
is obtained as discussed below. 

The three positions of the crank 
(0j, 0 2 and 0 3 ) may be obtained in 
the similar way as discussed in Art. 
25.10 and the corresponding three 
positions of the slider (s v s 2 and ,v ,j are 
obtained from the given condition as in 
equation (i). Now the dimensions a, b 
and c may be determined as discussed 
below : 



Fig. 25.20. Slider crank mechanism. 


In a right angled triangle BC'C , 

BC = b\ BC' = asin0 — c , and CC' = s — a cos0 


b 2 = (a sin 0 - c) 2 + 0 - a cos 0) 2 

9.9 9 . 222 

= a sin 0 + c —2acsin0 + 5 +a cos 0 — 2sacos0 


= a 2 +c 2 — 2czcsin0 — s 2 — 2ascos0 
or 2ascos0 + 2acsin0 + £> 2 -a 2 -c 2 = s 2 

/qscos0 + k 2 sin0-k 3 = s 2 
k\=2ci ; k 2 =2ac and k 3 = a 2 -b 2 +c 2 


... (tt) 


where 


... (hi) 
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For the three different positions of the mechanism i.e. for ( 61,62 63 ) and (,v p s 2 , s 3 ), the 
equation (ii) may be written as 

/q .v, cos 6 , +k 2 sin 6 , — k 3 = .v , 2 ... (iv) 

k x s 2 cos 6 2 + k 2 sin 6 2 - k 3 - s\ ... (v) 

k x s 3 cos 63 + k 2 sin 6 3 -k 3 -s 3 ... (vi) 

The equations (iv), (v) and (vi) are three simultaneous equations and may be solved for 
three unknowns k v k 2 and k?. Knowing the values of k v k 2 and k?, the lengths a, b and c may be 
obtained from equations (Hi). 

Example 25.9. Synthesize a slider crank mechanism so that the displacement of the slider 
is proportional to the square of the crank rotation in the interval 45° < 0 < 135° ■ Use three preci¬ 
sion points with Chebyshev’s spacing. 


Solution : Given. 6 S =45° ; 8 F =135° 

First of all, let us find the three precision points (i.e. x v x 2 and xf. We know that 


Xj = ^(x s +x F )-^-(x F -x s )cos 


K(2j-D 


2 n 


; where j = 1,2 and 3 


Assuming the starting displacement of the slider (s s ) = 160 mm and final displacement of 
the slider (s F ) = 30 mm. It may be noted that for the crank rotating in anticlockwise direction, the 
final displacement will be less than the starting displacement. 


x, = -(100 + 30)--(30-100)cos 
2 2 


71(2x1-1) 

2x3 


= 95.3 mm. 


...(•.• x s = s s ; jc f = s F and n = 3) 


x 2 =^(100 + 30)-^(30-100)cos 

and *3 = - (100 + 30) - - (30 -100) cos 

2 2 

The corresponding three values of 0 are given by 

0 — 0 

0 . = 0g H— - — (Xj — x s ) ; j = 1, 2, and 3 

Xp-Xg 

6 , =45+ 135 ~ 45 (95.3-100) = 51.04° 

1 30-100 


0. =45 + ——- (65-100) = 90° 

2 30-100 

and 0. =45 + 135 ~ 45 (34.7-100) =128.96° 

3 30-100 

Since it is given that the displacement of the slider (s) 
is proportional to the square of the crank rotation (0 ), therefore, 
for the displacement from initial position (s s ) to s when crank 
rotates from initial position ( 0 S ) to 0 , we have 


7t(2x2 —1) 
2x3 

7t(2x3-l) 

2x3 


= 65 mm 

= 34.7 mm 



A belt-conveyor that can trans¬ 
port small components. 

Note : This picture is given as 
additional information and is not a 
direct example of the current chapter. 
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s — ,? s = C (0 — 0 S ) 2 ... ( 9 is expressed in degrees) 

s-s s _ 30-100 _ -7 

”(0-0 s ) 2 “(135-45) 2 “8l0 - (Takin S * = *f ^ and 0 = 0 f) 

Now the three positions for the slider displacement ( s ) corresponding to the three positions 
of the crank angle (0 ) are given by 

s, =So+C(0, -0A 2 =100-— (51.4-45) 2 =99.7 mm 
1 s i s gl0 

s 2 =j s +C(0,-0c) 2 =100- — (90-45) 2 =82.5 mm 

2 S 2 S gl0 

s, = So +C(0, -0c) 2 = 100—^ (128.96 — 45 ) 2 = 39.08 mm 
3 s 3 s 810 

Now the three equations relating the (Oj,^), (0 2 ,s 2 ) and (631 5 3 ) are written as 
*! x 99.7 cos 51.04° + k 2 sin 51.04° - * 3 = {99.1) 2 


or 

62.7 kj + 0.7776 k 2 -k 2 = 9940 

... (i) 


Similarly, k i xS2.5cos90° + k 2 sm90°-k 3 =(82.5) 2 


or 

k 2 -k 3 = 6806 

... (ii) 

and 

k { x39.08cos 128.96° + k 2 sin 128.96° -k 3 = (39.08) 2 


or 

-24.57^+0.776jk 2 -jfe 3 =1527 

... (iii) 


The equations (i), ( 11 ) and (iii) are three simultaneous equations in three unknowns k v k 2 
and k y On solving, we get 

k x = 96.4 ; k 2 = 13 084 ; and k 3 = 6278 
We know that k\ = 2a, or a = k l / 2 = 96.4 / 2 = 48.2 mm Ans. 

k 2 = 2a.c or c - k 2 / 2a - 13 084 / 2 x 48.2 = 135.7 mm Ans. 

and k 3 =a 2 -b 2 +c 2 

b 2 = a 2 +c 2 -k 3 = (48.2) 2 + (135.7) 2 -6278 = 14 460 
or b = 120.2 mm Ans. 

EXERCISES 

1. In a four bar mechanism PQRS, the link PS is fixed. The length of the links are : PQ = 62.5 mm ; 
QR = 175 mm ; RS = 112.5 mm and PS = 200 mm. The crank PQ rotates at 10 rad/s clockwise. 
Find the angular velocity and angular acceleration of the links QR and RS for the values of angle 
QPS at an interval of 60°. 

2. In a slider crank mechanism, the crank AB = 100 mm and the connecting rod BC = 300 mm. When 
the crank is at 120° from the inner dead centre, the crank shaft has a speed of 75 rad/s and an 
angular acceleration of 1200 rad/s 2 both clockwise. Find at an interval of 60° 1. the linear velocity 
and acceleration of the slider, and 2. the angular velocity and angular acceleration of the rod, when 
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(a) the line of stroke of the slider is offset by 30 mm, and 

(b) the line of stroke of the slider is along the axis of rotation of the crank. 

3. A mechanism is to be designed to generate the function 

y = x ° 8 

for the range 1 < x < 3 , using three precision points. Find the three values of x and y. 

[Ans. 1.134, 2, 2.866 ; 1.106, 1.741, 2.322] 

4. Determine the three precision positions of input and output angles for a mechanism to generate a 
function 

y = x LS 

when x varies from 1 to 5, using Chebyshev’s spacing. Assume that the initial values for the input 
and output crank are 30° and 90° respectively and the difference between the final and initial 
values for the input and output cranks are each equal to 90°. 

[Ans. 36°, 75°, 94.48°; 91.22°, 144.57°, 181.22°] 

5. Synthesize a four bar linkage using Freudenstein’s equation to generate the function y = x 1 ' 8 for 
the interval 1 < x < 5. The input crank is to start from 0§ = 30° and is to have a range of 90°. The 
output follower is to start at (|><j = 0° and is to have a range of 90°. Take three accuracy points at x 
= 1,3 and 5. 

6. A four bar function generator is used to generate the function y = 1/x for 1 < x < 3 between the 
input angle of a crank and the angle the follower makes with the frame. Find the three precision 
points from Chebyshev’s spacing if the initial values of input angle (i.e. crank angle) and output 
angle (i.e. follower angle) are 30° and 200° respectively. The difference between the final and 
initial values of the crank and follower angles are each equal to 90°. 

7. Synthesize a four bar linkage that will generate a function y = x 1 2 for the range 1 < x < 5. Take 

three precision points : 9 S = 30°; <t>s = 60° and A0 = A<j> = 90°, where 9 S and (f> s represent respec¬ 
tively the initial angular positions of the input and output crank; A9 and A0 are respectively the 
ranges of the angular movements of the input and output crank. 

8 . Synthesize a four bar mechanism to generate the function y - log x, where x varies between 1 and 
10. Use three accuracy points with Chebyshev’s spacing. Assume 0 S = 45°; 0 F = 105°; (]) s = 135° 
and (|> F = 225°. Take the length of the smallest link equal to 50 mm. 

9. Synthesize a four bar mechanism to move the rod AB as shown in Fig. 25.21, through the positions 
1, 2 and 3. The end points A and B are used as moving pivot points. 



Fig. 25.21 



Fig. 25.22 
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10. Design a four bar mechanism to guide the door in and out with little rotation until it clears the 
surrounding structure, after which it swings fully open to one side. The three positions of such a 
door under going this type of motion is shown in Fig. 25.22. The points A and B are used as 
moving pivots that guides the body through the three positions. 


DO YOU KNOW ? 


1. Explain Freudenstein’s method of three point synthesis of mechanisms. 

2. Derive the expressions for displacement, velocity and acceleration of a four bar mechanism. 

3. What do you understand by coupler curves ? Describe the method of obtaining the co-ordinates of 
a coupler point in a slider crank mechanism. 

4. Explain synthesis of mechanism with examples. What do you understand by 

(a) Type synthesis ; ( b ) Number synthesis ; and (c) dimensional synthesis. 

5. Describe the classifications of synthesis problem. 

6. Write an expression for determining the precision points. 

7. Discuss the method of determining the angles for input and output link in a four bar mechanism for 
function generation. 

8. Describe the method of designing a four bar mechanism as a function generation. 


OBJ ECHVETYPE QUESTIONS 

1. The analysis of mechanism deals with 

(a) the determination of input and output angles of a mechanism 

( b ) the determination of dimensions of the links in a mechanism 

(c) the determination of displacement, velocity and acceleration of the links in a mechanism 

(d) none of the above 

2. The synthesis of mechanism deals with 

(a) the determination of input and output angles of a mechanism 

( b ) the determination of dimensions of the links in a mechanism 

(c) the determination of displacement, velocity and acceleration of the links in a mechanism 

( d ) none of the above 


The three precision points in the range 1 < x < 3 are 

(a) 1.1, 2, 2.6 ( b ) 1.6, 2.5, 2.95 

(c) 1.134, 2, 2.866 ( d ) 1.341, 2,2.686 

For a four bar mechanism, as shown in Fig. 25.23 the Freudenstein’s equation is 

(a) kjcosB + ^costjt + A^ = cos (0 — <()) 

( b ) kjcos0-k 2 cos(]) + k 3 = cos(0-(])) 

(c) k!cos(|) + k 2 cos 0 + k3 =cos(0 —(|>) 

(, d ) kjeostIt-kjeosB + Ay =cos(0-(])) 


d d j _a — b +c +d 

where h~~ ; * 2 = — ; *3 = 


lac 

ANSWERS 



1. (C) 


2. (Z>) 


3. (c) 


4. {d) 



To FIRST 





